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INTRODUCTION 

The  term  fluidlzation  Is  used  to  describe  a  suspension  of 
solid  particles  in  a  fluid  flowing  through  a  container.   This 
combination  of  particles  and  fluid  constitute  a  fluidlzed  bed. 
Within  this  general  concept  of  fluidlzation  there  are  many  spe- 
cial cases  which  arise.   First,  the  fluid  employed  may  either  be 
liquid  or  gaseous.   Each  of  these  present  individual  problems  of 
analysis.   The  solid  being  fluidlzed  might  either  be  more  dense 
or  less  dense  than  the  fluidizing  fluid.   This  first  case  would 
give  rise  to  the  familiar  situation  of  the  fluidizing  medium, 
flowing  upward  through  the  bed.   The  latter  case  would  envolve  a 
downward  flow  of  fluid  and  has  never  been  treated  to  the  author's 
knowledge.   Two  distinct  mechanisms  of  fluidlzation  have  been  de- 
fined.  First,  there  is  particulate  fluidlzation  in  which  each 
solid  particle  is  considered  to  act  as  an  individual  unit  within 
the  bed.   Then  there  is  aggregative  fluidlzation  in  which  the 
particles  are  not  considered  as  individual  units  but  as  part  of 
groups  or  aggregates.   The  distinction  between  aggregative  and 
particulate  fluidlzation  was  first  observed  by  Wilhelm  and  Kwauk 
(1).   Particulate  fluidlzation  is  characterized  by  the  existence 
of  a  homogeneous  particle  distribution,  while  in  aggregative 
fluidlzation,  two  phases,  dilute  and  dense,  are  seen  to  exist. 
Finally,  the  geometry  of  the  container  may  be  chosen  at  will. 
Most  of  the  work  to  date  has  been  carried  out  in  columns  with 
circular  cross  sections,  but  there  have  been  a  few  studies  made 
using  columns  with  square  cross  sections. 


This  work  is  concerned  with  the  fluidizatlon  of  soherical, 
solid  particles  by  a  liquid  of  less  density  in  a  column  of  circu- 
lar cross  section.   The  very  nature  of  the  experimental  and  ana- 
lytical work  envolved  imply  that  the  fluidizatlon  is  to  be  partic- 
ulate  or  very  nearly  particulate. 

The  use  of  fluidizatlon  in  industry  is  very  wide  and  is  in- 
creasing in  scone  (2).   A  few  of  the  more  important  applications 
are  summarized  as  follows: 

1)  Chemical  reactors. 

a)  Where  the  fluid  stream  is  composed  of  the  reacting 
component  and  the  solid  ^articles  catalyze  the 
reaction. 

b)  Where  reaction  takes  place  between  the  fluid  8nd 
solid. 

2)  Mass  transfer  operations. 

a)  Leaching. 

b)  Adsorotion. 

c)  Drying. 

3)  Evaporation. 

a)   Solid  particles  can  be  fluldlzed  in  a  boiler  to  pre- 
vent scale  formation  on  heat  transfer  surfaces, 
M   Nuclear  reactors. 

a)  Fluid  ^ctins-  as  moderator  and  ^articles  as  fuel  ele- 
ments. 

b)  Fluid  acting  as  fuel  element  and  narticles  as  mod- 
erator. 


5)  Elutriation  or  separation  of  particles  by  their  different 
fluldizing  characteristics. 

6)  Calcination  and  heat  treatment  processes. 
Generally  the  fluidlzation  technicue  will  yield  relatively 

high  heat  and  mass  transfer  coefficients  and  will  provide  a  high 
degree  of  mixing  between  fluid  and  particles. 

Process  Control  Considerations 

In  many  of  the  applications  summarized  the  maintenance  of  a 
given  oorosity  or  bed  height  is  either  desirable  or  essential. 
For  instance,  the  neutron  density  of  a  fluidized  nuclear  reactor 
will  be  a  function  of  porosity,  and  failure  to  maintain  a  given 
rorosity  within  certain  limits  will  result  in  a  non-critical  con- 
dition.  Also  in  any  application  of  fluidlzation  there  will  be  a 
porosity  determined  by  design  procedures  or  by  experimental  work 
at  which  the  orocess  is  to  operate.   Therefore,  the  problem  of 
process  control  is  becoming  increasingly  apparent  in  the  area  of 
fluidlzation  as  in  all  phases  of  the  chemical  process  industry. 
In  tne  past,  a  large  portion  of  process  control  has  been  accom- 
plished by  human  operators  directly  observing  the  process  and  mak- 
ing necessary  corrections.   At  the  oresent  time  the  need  for  pro- 
ducts of  higher  quality  and  the  high  speed  response  of  equipment 
have  rendered  the  human  approach  to  control  Impractical.   Finan- 
cial considerations  also  dictate  a  need  for  automatic  control  of 
processes.   Ceaglske  (3)  defines  automatic  control  as  follows: 

"A  control  system  consists  of  a  series  of  units  combined 
to  produce  a  desired  effect  with  little  or  no  human  supervision. 


The  basic  parts  of  the  oynt3m  are  the  process,  the  measuring  ele- 
ment, the  error-detecting  mechanise,  thn  controller,  «nd  the  final 
control  element." 

The  measuring  element  continuously  mearures  the  value  of  some  pro- 
cess  variable  and  sends  a  signal  which  is  proportional  to  the 
measured  variable  to  the  error-detecting  mechanism.   The  error- 
detecting  mechanism  then  compares  the  value  of  the  measured  vari- 
able with  the  desired  value  and  indicates  this  difference,  or 
error,  to  the  controller.   The  controller  takes  corrective  action 
based  on  the  error  signal  and  indicates  to  the  final  control  ele- 
ment the  necessary  change  to  be  effected  in  the  input  variable. 
Therefore,  the  system  consists  of  a  do  ed  circuit  in  which  the 
lnnut  ftt  any  time  is  dependent  ur>on  the  output  at  a  time  Ju?t 
preceding. 

The  nature  of  the  corrective  action  taken  depends  uoon  the 
type  of  controller  emoloyed.   A  proportional  controller  is  one 
which  dictates  a  corrective  action  pro  ortlonal  to  the  error  de- 
tected.  The  Integral  control"!  er  BOl     m   the  integral  of  the 
error  with  time.   Likewise,  ther~  are  rat~  dApendent  and  two-po- 
sition controller  as  well  «is  combinations  of  these  various  modes. 

Each  of  the  elements  of  a  control  tf     has  its  own  response 
characteristics,  and  the  nerformsnce  of  the  system  as  a  wh->la  Is 
dependent  upon  a  combination  o*"  the  various  individual  chpmcter- 
lstlcs.   It  is  therefore  necessary  for  the  system  designer  to 
have  at  his  disposal  Information  on   the  dynamic  characteristics 
of  each  component  of  the  system.   By  for  the  most  important  com- 
ponent is  the  process  itself  si  nee  the  design  of  all  the  other 
component?  is  dependent  on  the  dynamic  behavior  of  this  •lament. 


The  disturbances  a  process  may  experience  are  usually  random 
in  nature,  and  it  is  very  difficult  to  determine  the  resoonse  of 
the  system  to  random  inputs.   Hov:evert  it  ie  possible  to  describe 
the  response  of  a  system  to  any  input  as  a  combination  of  its  re- 
sponses to  cert.ln  regular  inputs.   In  particular  the  transient 
response  to  a  step  input  and  the  steady-state  response  to  a 
sinusoidal  input  are  useful. 

In  analyzing  the  dynamic  characteristics  of  a  system  one  of 
two  methods  is  generally  used.   The  first  involves  the  deriva- 
tion of  a  differential  equation  describing  the  process  and  solv- 
ing it  fo^  the  desired  inputs  to  obtain  response  characteristics. 
This  method  is  quite  convenient  when  the  differential  equation  ad- 
mits an  analytical  solution  in  closed  form.   However,  when  the 
equation  requires  solution  by  series  or  numerical  methods,  the 
solution  obtained  is  not  nearly  as  useful  for  control  system  de- 
sign.  A  relatively  new  and  very  pov.-arful  raeth-d  of  system  analy- 
sis Involves  the  use  of  analog  computers.   Use  of  this  method  re- 
quires only  that  a  differential  equation  describing  the  process 
be  available.   Electrical  circuits  in  the  luaalftj  are  set  up  ac- 
cording to  the  equation,  and  the  analog  computer  than  simulates 
the  process.   Any  input  can  be  fed  Into  the  computer  and  the  re- 
sponse noted.   This  information  is  then  used  in  the  design  of  the 
actual  control  system.   This  method  of  analysis  Is  particularly 
useful  in  the  analysis  of  non-linear  systems. 

The  fluldlzed  bed  presents  a  typical  non-linear  problem.   The 
non-linear  element  existing  in  fluldlzed  bed  analysis  can  be  con- 
sidered roughly  analogous  to  the  non-linearities  occurlng  in  heat 
and  mass  transfer.   Just  as  the  heat  and  mass  transfer  coeffi- 
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clents  are  dependent  upon  the  fluid  flow  conditions,  the  relation- 
ship between  fluldizing  velocity  and  porosity  is  also  dependent  on 
flow  conditions.   Although  a  uite  a  bit  of  work  has  been  done  on 
the  analysis  of  the  dynamic  behavior  of  heat  and  mass  transfer 
equipment,  little  has  been  done  on  the  fluid i zed  bed  expansion. 

Purpose 

The  major  purpose  of  this  work  was  to  develop  a  linearized 
model  to  describe  the  dynamic  behavior  of  the  liquid-solid  fluld- 
lzed  bed  and  to  verify  that  model  with  experimental  results.   It 
was  hoped  that  a  sufficiently  accurate  representation  could  be  ob- 
tained to  be  useful  in  the  design  of  control  systems  for  fluidl- 
zation  processes. 

EXPERIMENTAL 

Apparatus  and  Equipment 

Photograohes  of  the  equipment  are  shown  in  Pigs.  1  and  2  and 
a  schematic  diagram  is  shown  in  Fig.  3» 

The  column  used  was  a  ?ii~inch  high,  1.982-inch  diameter  glass 
column.   The  column  was  seated  on  a  niece  of  J/k^i.nch   Lucite  with 
a  2-inch  diameter  hole  through  the  center.   A  circular  Piece  of 
U.  S.  Standard  200  mesh  screen  was  glued  in  the  center  of  the 
Lucite  slab.   A  piece  of  U.  S.  Standard  30  mesh  screen  was  glued 
directly  under  the  200  mesh  in  order  to  provide  support.   This 
screen  assembly  formed  the  bottom  boundary  for  the  bed  and  also 
provided  a  degree  of  mixing  or  distribution  for  the  fluldizing 
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stream.   The  column  and  its  accessories  sat  directly  over  a  calm- 
ing section.  The  bottom  portion  of  the  calming  section  was  six 
inches  high  and  had  an  inside  diameter  of  four  inches.   The  top 
portion  which  was  connected  to  the  Lucite  screen  supporter  was 
3  1/2  inches  high  with  a  two  inch  inside  diameter.   The  calming 
section  was  packed  with  1/4-inch  Berl  saddles  to  such  a  height 
that  the  top  of  the  packing  was  1/8-  to  l/^-inch  below  the  sup- 
porting screens.   This  packed  region  was  used  for  the  purpose  of 
providing  a  more  or  less  flat  velocity  profile  at  the  entrance  of 
the  column.   The  fluid  entered  the  calming  section  through  a  nip* 
welded  in  its  side  at  the  bottom.   At  the  top  of  the  column  anoth- 
er 200  mesh  screen  was  held  in  place  between  two  rubber  gaskets. 
This  prevented  any  of  the  smaller  particles  from  being  carried  out 
of  the  column  and  through  the  system.   Immediately  above  this 
screen  a  short  section  of  two  inch  pipe  was  placed.   This  oipe 
was  equipped  with  a  temperature  well,  thus  providing  a  means  of 
reading  the  temperatures  of  the  liquid  Just  as  it  left  the  bed. 
The  water  supply  was  a  50  gallon  constant-head  tank.   Water  which 
had  circulated  through  the  system  was  returned  from  the  top  of  the 
column  back  to  the  supply  tank.   The  only  problem  encountered  in 
this  scheme  was  the  slight  rise  in  temperature  the  water  experi- 
enced upon  flowing  through  the  system.   Since  the  experimental 
work  was  to  be  conducted  under  isothermal  conditions,  it  was 
necessary  to  continually  add  a  small  flow  of  tap  water  to  the 
supply  tank  and  let  the  excess  overflow  through  a  pipe.   This 
procedure  worked  out  very  well  and  it  was  possible  to  maintain 
a  constant  temperature  within  ±  lop.   The  water  was  forced  through 
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the  system  by  means  of  an  Eastern,  Model  P,  pump.   The  flow  through 
the  system  was  measured  by  means  of  two  rotameters  in  parallel. 
Only  one  of  these  rotameters  was  utilized  for  any  one  experiment. 
The  choice  deoended  on  the  desired  range  of  flow  rates.   Calibra- 
tion curves  for  the  two  rotameters  are  shown  in  Figures  4- (a)  and 

4-(b). 

When  it  was  desired  to  produce  a  step  or  pulse  input  in  flow 
rate,  the  flow  was  directed  through  a  circuit  consisting  of  two 
solenoid  valves  in  parallel.   Each  of  these  valves  was  in  series 
with  a  globe  valve  which  could  be  used  to  regulate  the  flow 
through  that  particular  line  when  it  was  open.   The  solenoids 
used  were  Sporlan,  Type  12  P  and  were  3/8  inch,  two-way  valves. 
They  were  connected  through  one  electrical  switch  in  such  a  man- 
ner that  one  would  be  open  and  the  other  closed  on  each  switch  po- 
sition. 

Sine-wave  generator:   The  apparatus  used  for  oroducing  a  sine 
or  cosine  input  in  fluldizing  velocity  co- sisted  of  five  main  u- 
nits: 

(1)  Power  for  the  apparatus  was  supplied  by  a  General  Elec- 
tric, 1/2  h.p*,  Model  5KC63AB66B,  A.C.,  motor  operating 
at  1725  R.P.M. 

(2)  The  outout  from  the  motor  drove  a  Vickers  3/4  h.p.  hy- 
draulic transmission,  Model  No.  TR3-HR13-F18-13,  with  a 
maximum  input  and  output  speed  of  1800  R.P.M,   This  was 
a  variable  speed  transmission  giving  an  output  ranging 
from  0  to  100^  of  the  input  speed. 

(3)  The  output  from  the  Vickers  transmission  was  further  re- 
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duced  by  use  of  a  200-1  reduction  gear  assembly.   This 
was  necessary  due  to  erratic  performance  of  the  hydrau- 
lic transmission  at  the  very  low  frequencies  desired. 

(4)  The  conversion  from  the  constant  speed  rotational  output 
from  the  gear  assembly  to  longitudinal  motion  with  a  sine 
varying  velocity  was  accomplished  usir>g  the  apparatus 
described  in  Pig.  5.   The  lever  arm  shown  was  equipped 
with  drilled  and  taoped  holes  at  1/k   inch  intervals,  thus 
providing  lengths  from  1-1/2  Inches  to  5-3 /k   inches, 

(5)  The  final  unit  was  a  Model  2k   Alkon  Cylinder  with  a  12 
inch  maximum  stroke.   The  diameter  of  this  cylinder  was 
1-3/4  inches. 

A  photograph  of  the  sine-wave  generator  is  shown  in  Figure  6. 

Particles 

Five  different  groups  of  particles  were  employed  in  this  ex- 
periment.  The  particle  size  range  within  each  group  and  the  par- 
ticle density  were  as  follows: 

(I)   These  ^articles  were  Propper  Borosilicate  size  3  beads, 
6  to  7  U.S. -mesh.  Particle  diameters  ranged  from  2.83  nim. 
to  3«36  mm;  the  average  diameter  was  2.8?  mm.   As  is  evi- 
dent by  observing  the  average  diameter  in  comparison  with 
the  two  extremes,  the  size  distribution  was  not  normal. 
In  ^"act,  25  of  these  particles  were  measured  with  a  mi- 
crometer giving  an  average  of  2.87  mm.  and  a  variance  of 
only  2.37  x  10"^mm2.   The  density  of  the  particles  was 
2.5  g./cc. 
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Let        Vs  =  Velocity     of    shaft     to     cylinder 
VL  =  Tangential     velocity    of    lever    arm 
&  =  Angle    between    Vs    and     VL 
6  =  Angle    between     lever    arm     and    shaft     to    cylinder 
Then,  Vs  =  VLCosCX 

^  =  90-  0< 
Vs  =  VLSin£ 

Since    the    cross- sectional    area    of     the    cylinder    is    constant 
the    flow    rate     produced     by    a     constant    rotational 
velocity,    VL;    is 

V  =AcyVs  =AcyVL  Sin^ 
Where         V  =  Volumetric     flow     rate    (ml/sec) 

Acy=  Cross    sectional     area     of     cylinder     (cm2) 

Fig.  5.    Assembly     for     sine-wave      production. 
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(II)  These  particles  were  Nalcite  ion  exchange  resin  manu- 
factured by  the  Nalco  Chemical  Company.   The  particles 
were  screened  and  the  25  to  30  U.S. -mesh  fraction  used, 
giving  a  particle  diameter  range  of  0.59  mm.  to  0.71  mm. 
Since  this  material  absorbs  a  considerable  amount  of 
water,  all  screening  and  physical  property  estimation 
were  conducted  with  the  particles  saturated  with  water. 
The  density  of  the  ^articles  was  1.08  g./cc. 

(III)  These  particles  as  well  as  those  in  particle  groups  IV 
and  V  were  glass  beads  manufactured  by  the  Minnesota 
Mining  and  Manufacturing  Company.   This  group  consisted 
of  the  ^0  to  fcj  U.S. -mesh  fraction  giving  a  particle 
diameter  range  of  0.35  mm.  to  0.42  mm.   The  density  of 
these  particles  war  2.5  g./cc. 

(IV)  These  particles  were  identical  to  those  of  group  III 
with  the  exception  that  they  were  the  60  to  70  U.S. -mesh 
fraction  giving  a  particle  diameter  range  of  0.21  mm. 

to  0.25  mm. 

(V)  These  particles  were  identical  to  those  of  groups  III 
and  IV  with  the  exception  that  they  were  the  70  to  80 
U.S. -mesh  fraction  giving  a  particle  diameter  range  of 
0.177  mm.  to  0.21  mm. 

Henceforth,  these  particle  croups  will  be  specified  by  their  U.  S. 
mesh  range. 

The  close  cuts  in  particle  size  range  were  considered  necessary 
in  order  to  minimise  the  effect  of  si.^e  distribution.   Another  as- 
pect which  must  be  considered  is  the  sphericity  of  the  particles. 
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In  all  cases   the  oarticles  are  c onsldered   spherical  and    the  fol- 
lowing evidence   is   given   for  this  assumption.      In  group  I,    the 
particles  were  large   enough  to  permit  measurement  by  micrometer. 
Therefore   several   of  these  particles  were  measured  along  different 
axes  to  obtain  an  indication  of  their  symmetry.      The   following 
data  are  representative  of  these  measurements: 


Table  1,      Measurement  of  particles   in  Group   I 
for   symmetry. 


Particle   No.  Diameter    (inches) 


1  0.1132 

"  0.1129 

"  0.1133 

■  0.1137 


2  0.1113 
•  0.1129 

0.1110 

"  0.1128 

3  0.1163 

1  0.13.68 

"  0.11^8 

"  0.1158 


The  oarticles   In   the   other  groups  were  too   small  to  allow  in- 
dividual measurement.      In  this   case  an   indication  was   obtained  by 
comparing  the  values   of  terminal   falling  velocity   observed  exper- 
imentally with  those  calculated  assuming  spherical  particles. 
These  figures  are  presented   in  Table  2. 
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Table  2.   Terminal  velocities  of  particles. 


Particle  Group  No.  :  Ug  (calculated)  :  Ug  (experimental) 
: (cm./rec.)  : (cm. /gee.  ) 

II  l.oo  1.15 

III  4.82  5.06 

IV  2.68  2.83 

V  2.12  2.35 


It  can  be  seen  that  in  all  cases  there  is  fair  agreement,  and 
while  this  analysis  does  not  provide  a  quantitative  prediction  of 
sphericity,  it  does  indicate  that  the  particles  have  a  sphericity 
close  to  1.0. 

Procedure 

General:   The  desired  weight  of  particles  of  a  particular 
group  were  placed  in  the  column.   The  temperature  of  the  fluidiz- 
ing  water  was  fixed  and  maintained  at  some  constant  value.   Steady 
state  fluidization  data  were  then  obtained  by  setting  the  flow 
rate  at  a  constant  value,  allowing  a  sufficient  amount  of  time  for 
the  bed  to  reach  equilibrium,  and  reading  the  rotameter  and  bed 
height.   All  bed  height  readings  were  made  using  a  scale  calibra- 
ted in  centimeters  which  was  pasted  on  the  side  of  the  column. 
This  was  repeated  for  various  bed  expansions  and  the  results  re- 
corded. 

Step  Chan ge  Data :   The  bed  was  allowed  to  come  to  and  remain 
at  some  steady-state  height.   A  switch  was  thrown  which  caused 
both  solenoid  valves  to  operate  and  effect  a  change  in  flow  from 
one  line  to  the  other,  giving  a  step  change  in  flow  rate  to  the 
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bed.   The  instantaneous  bed  height  was  read  at  regular  time  in- 
tervals during  the  transition  from  the  initial  steady-state 
height  to  the  final  steady-state  height. 

Pulse  Data:   The  response  to  a  pulse  composed  of  two  step  in- 
puts was  obtained  in  much  the  same  manner  as  the  response  to  a 
step  change.   The  only  difference  in  procedure  was  to  effect  a 
second  step  change  in  flow  back  to  the  original  flow  rate  at  some 
time  before  the  bed  height  reached  the  steady-state  value  corres- 
ponding to  the  intermediate  flow  rate. 

Frequency  Response  Da bat   A  steady  state  flow  rate  was  first 
set  and  the  bed  allowed  to  come  to  equilibrium.   The  sine  wave 
generator  was  adjusted  for  lever  arm  length  and  rotational  veloc- 
ity in  order  to  obtain  the  desired  input  amplitude  and  frequency. 
The  sine  input  was  then  superimposed  on  the  steady-state  flow  rate, 
After  waiting  a  few  cycles  to  insure  that  the  transient  effects 
had  sufficiently  subsided,  the  maximum  and  minimum  in  bed  height 
were  read  and  recorded.   Phr-se  lag  data  were  obtained  by  observ- 
ing and  recording  the  time  difference  between  the  reversal  in 
piston  direction  and  the  reversal  in  bed  direction.   This  was  not 
possible  when  operating  at  high  frequencies  where  the  phase  lag 
was  close  tc  90°  and  the  time  interval  between  reversals  was  very 
small. 

Effect  of  Load  Changes  on  Steady  State 
Fluidizlng  Velocity 

One  question  which  might  arise  pertaining  to  the  experimental 
procedure  is  the  effect  of  the  transient  flow  producing  elements 
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of  the  system  on  the  performance  of  the  supply  pump  and  therefore 
on  the  steady-state  flow  delivered.   Figure  3  snov.'s  that  these 
elements  are  directly  connected  to  the  steady-state  circuit. 

In  particulate  fluldlzation  the  pressure  drop  through  the  bed 
is  very  nearly  independent  of  flow  rate  or  bed  heigth.   Therefore, 
for  the  step  input,  the  only  change  in  load  is  due  to  increased 
friction  losses  in  the  lines.   It  is  assumed  that  this  effect  is 
almost  instantaneous  so  that  little  distortion  of  the  ir.put  occurs. 
For  the  sine  wave  input  the  situation  is  quite  different.   The 
sine  generator  is  imposing  a  continually  changing  load  on  the  sys- 
tem and  it  might  be  expected  that  some  of  its  effect  would  be 
transmitted  back  through  the  system  causing  the  steady-state  flow 
to  vary.   It  would,  however,  be  possible  to  detect  any  s uch  effect 
by  observing  the  rotameter.   No  oscillating  of  the  rotameter  was 
observed  during  any  of  the  experiments,  thus  confirming  that  a 
valid  steady-state  flow  was  being  obtained*   This  was  evidently 
due  to  the  fact  that  the  load,  changes  were  small  and  that  the 
pump  was  operating  on  a  fairly  level  portion  of  Its'  lo^d-outout 
curve.   At  sinusoidal  inputs  of  greater  amplitude  than  were 
necessary  in  this  experiment,  oscillations  were  observed  in  the 
steady-state  flow  and  this  is  definitely  a  problem  in  any  experi- 
mentation of  this  nature. 

DEVELOPMENT  OF  LINEAHI^ED  MATHEMATICAL  MODEL 

As  was  stated  in  the  Introduction,  the  major  analytical  tool 
used  in  this  work  was  a  linearized  description  of  the  basically 
non-linear  process. 
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Richardson  and    "aki    (4)    showed  that  for   steady-state  behavior 
the  porosity  can   be  related   to   the  fluidizlng  velocity  by 

4>ss   =      Ua€^  (1) 

where    ^Ss     ls   the   steady-state  superficial    liauld   velocity,   Ug   is 
the  terminal  falling   velocity   of   the  particle,  £  ss     is   the  steady- 
state  porosity,    and  n   is  a   constant   characteristic   of   the  parti- 
cle-llquid    system. 

In   considering  the  response  of  a   liouid-solid   fluidized  bed 
to  a    ste-n-innut    in   fluidizlng   velocity,    Slls,    et.    al.    ( c )    extend- 
ed  Equation    (l)    to   the  unsteady-state  case  and   obtained  the   fol- 
lowing expression  for   the  rate  of  chance  of  the  bed  height: 

•  -'♦**  V-.  (?) 

where  h   is   the  instantaneous  bed  height,    t   is   the   time,  rl     is 
the  superficial   liquid  velocity  after  the  step   input,   and      C.      is 
the   instantaneous  porosity  at   the   top  of   the  bed. 

For  the   steady-state  case,    dh/dt    =      0,    4>  i   =     <^g8    ,    and 
Equation    (?)   reduces   to   Equation    (1). 

Equation   (2)   can   be  linearized  by   letting 

*1      -    4>ss+     *+  0) 

h  =    hsg     ■+       /Jh  (*) 

eh     -   e8S  +  ^e  (5) 


Mow   £  **  can  fee   expanded  In   a    m^;lor  eerier   about    £  ss  tc  give 

n  n  r.-l  M  ,  z  v 

€  _      -=.    e         +     n   €        4€  +  -    -  -  (6) 

h  ss  ps 

Substitution  of  Equations  (3),  (*0  and  (6)  into  Equation  (2)  gives 

(7) 


Since  He   rmount   of   eclids   in  thi    bed  rerr^rr    crnstart,    cn  1  n- 
rtanteneous  nieterlrl  balance  is  vritter   ai 

(1  -e)h    -     (l  -    ess)  hss  (8) 

where  G  is  the  instantaneous  porosity. 
Differentiation  of  Fcuatlon  (8)  fives 


iL  -   (1"  £BS)hSS 

dh  h2 


(9) 


for  small   changes   in  h,    Equation    (9)    becomes 


Ae      -    |*f  €U)   A  h  (10) 

hss 


Substituting  Sanation    (10)    into   Equation    (7)    <?ives 


d(hDO)   dMh)  ,  1-  e. 

'ss 


dt       dt  s   ss    s    &s   h 


but    dh00 

dt  '  ss       8   ^  SS 


2k 

rrd   consequently  Equation    (11) 

reduces 

to 

UA*1   + 

1 

An 

<d  <p 

(12) 

dt 

ss 

s       "-ss      v          ss' 

<f>    can  be 

arbitrarily 

written 

as 

1 

-     *cs4 

4> 

(13) 

hss 

»s  ■  ^S1 

<>-€..> 

Substitute 

.en  of   Ecuat5 

on    (13) 

into    Equation    (1?) 

■fl  6    a   1 

shange  of 

dependent 

variable  frc 

un  A  h   tc 

i  h  pi  ves 

* 

dh  + 
dt 

1 

h 
ss 

h 

At 

(14) 

hss 

uB  mef*   (1- 

•  e     ) 

ss 

Us  n  6JJ*   (1- 

e    ) 

ss 

If  we 

consider  the 

I  liquid- 

•  solid   fluidized   bed 

to  be 

I   first 

order   (time  constant^ 

systen, 

the  rate 

of  change 

of  height  can  be 

written  at 

dh            1 

H     -h 

dt            T 

m     h. 

T 

(15) 

where  T   is 

:   the   time  constant  for   the  system  and   h 

^   is    the  input 

function. 

For  a 

step  change 

in  fl  uidizinp  vel 

.ocity,    the 

input 

function, 

hj_,    is  a   constant,    A-, 

so   that 

Equation 

(15)    beconn 

3S 

25 

dt            T               T 

(16) 

Comparing  Equation   (14)  with  Equation 

(16)    si: 

ows 

that  the 

| 

time   constant   of  the  system   Is 

T      -                      h|| 

(17) 

1                       n  1 

U  n  €  -1  x   (1-  e      ) 
s          ss              se' 

By 

substituting  Equation    (1)    into   Equation    (17), 

the 

time   con- 

stant  can  also  be  written  as 

rp          SS          ^     SS 

(18) 

tSsn     d-£8S> 

Either  Equation    (17)    or  Equation    (18)   can 

be  used 

to 

calculat 

;e 

the 

theoretical    time  constant   of  the  llqu 

id-solid 

fluidized  bed. 

Eelvtieni  for  Standard 

Tr   uts 

The  solutions  of  Equation    (35)   or   (16 

)   for  variou 

s   input 

functions   are  prlven  below. 

Step   Innut:      Integration  of   Equation 

(15)   or 

(16) 

gives, 

for 

AL  = 

A  -  Step  size, 

(19) 

or 

h"hss     -=      1   -    B»l/V 
A 

(20) 

as 

a    solution  for   the   step   inout. 
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"'vise   Imut :  The  Input  considered  here  is  composed  of  two 
ste-o  inputs  separated  by  a  time  interval,  a.   The  second  step  is 
of  the  same  magnitude  as  the  first,  but  opposite  in  direction. 
The  solution  to  Equation  (15)  for  t  his  inout  is 

^L  ,   (i.e-t/T)  m   u(t_a)   e    T  (21) 

A 

(t-a)  >  0 
where  u  is  the  unit  itap  function: 

u(t-a)  =:   0   for  0  <  t  <  a 

=  1   for  a  <  t 

^amr>  In  nut:   The  resnome  to  a  linearly  varying  input,  hj_-Bt, 
is 

h-h 


ss  _.   BT  (e-t/T   t  _  1}  (??) 

T 


Sinusoidal    Input :      For  the    input  h.  -   h*   sincot,    the   solution 
to  Equation    (15)    is 

h~hss   _       Bin(  co  t-h°<   ) 

-     - rr  (23) 

h*  Vli-T2  co    2 

where  h     it   the  input   amplitude   in  terirc   of  the  height   variable, 
CO    is  the  ftngalftr  frequency,   and    of     is-   the  phase  lag. 

Prom  Equation    (?J)    it    is   seen   that   the  axplltude  ratio   is   given 
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by 


A#R#  __   output  amplitude  _  hV/lt  u3  ?T2 
Input  amplitude  h 


(24) 


/1+  co  2qi2 


Bid  the  ih'^e  lag, 

e><    ■=■   -arc  tan  oo  T  (25) 

Ml  ALIUS  OF  DATA 

As  stated  In  the  experimental  section,  the  experimental  data 
consisted  mainly  of  the  measurement  of  bed  height  as  a  function  of 
time  under  the  variable  flow  rate.   Flow  rate  and  properties  of 
the  liquid-solid  system  were  also  recorded.   Data  were  obtained 
on  the  system  response  to  three  inputs;  step  input,  pulse  input, 
and  sinusoidal  Input 4   For  analysis  the  data  obtained  for  each  of 
these  inputs  will  be  treated  separately. 

Response  to  Steo  Input 

The  data   taken   on  the   response  to  a   step-input   la  fluidizixig 
velocity  was    bo   be  used   i'.o   verify  the   linearized  model  developed 
in  the  previous    section.      Figure  7   shows   typical    step  change  re- 
sponse curves. 

Referring  to   Eouation    (19)    it   can  be   seen  that  a   plot   of 
time,    t,    versus   the  group   In        1   -  ■■  "  ,ss      should   yield  a    straight 
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line  of  slo-e,  1/T.   Therefore,  by  plotting  the  experimental  data 
expressed  in  this  form,  an  experimental  determination  of  the  time 
constant,  T,  w^s  obtained  for  each  experimental  run.   A  typical 
determination  of  the  experimental  time  constant  is  shown  in  Figure 
8.   It  should  be  noted  th^t  the  time  constant  is  also  given  by 


the  time  corresponding  to  In 


1  _  h-hss 


-  -  1,   This  is  immedi- 


A 

ately  apparant  from  the  form  of  Equation  (19).   Since  the  points 

on  Figure  8  fall  very  nearly  on  a  straight  line  for  the  first  35 
seconds,  there  is  very  little  question  as  to  how  the  straight  line 
is  to  be  drawn.   However,  for  various  reasons  to  be  discussed 
later,  the  experimental  point*  did  not  always  behave  in  this  man- 
ner.  When  curvature  ^t  the  origrin  was  observed,  the  straight 
line  wps  always  drawn  so  as  to  represent  the  tangent  to  the  curve 
at  the  origin.   This  procedure  was  adopted  because  of  the  manner 
in  which  the  linearized  model  was  developed.   When  h  —   hss, 

In   1  -  h-*nss   —  o,  so  th-?t  the  origin  represents  the  point 

A   - 
where  the  model  should  best  represent  the  response.   An  instance 

where  this  procedure  was  necessary  is  shewn  in  Figure  9.   The 
rapid  deviation  of  this  curve  from  the  straight  line  is  caused  by 
the  fact  that  the  curve  represents  a  response  to  a  step-down  in 
flow  rate  and  this  case  is  not  aoproximated  nearly  as  well  as  the 
response  to  a  step-up  by  the  linearized  model.   It  is  seen,  how- 
ever, that  the  approximation  is  accurate  for  a  short  period  of 
time. 

In  the  development  of  the  linearized  model  the  theoretical 
expression  for  the  time  constant  was   found  to  be 
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h  ss 


nUs  €£•*  (l.e J 


The  values  of  hSB  and  <E  sg  are  fixed  by  the  experimental  operating 
conditions  and  values  of  Us  were  determined  experimentally  and 
calculated  as  discussed  in  the  experimental  section.   The  equations 


used  to  calculate  Ug  were 


Stokes  Law:     Us  -  ?qDP  (  ^  8"  ^  L}    0.000l<  NHe<  2.0     (26) 

18 


intermediate  law:  U8  =  2^2  ^7\1'1^  '  H-  ^L>°-71 

^L0.29       y    0.*3 

2*°<NRe<-500  (27) 

where  gc  is  the  gravitational  constant,  /  s  is  the  density  of  the 
particles, /^l  is  the  density  of  the  liquid, >/  is  the  viscosity 
of  the  liquid,  and  %e  is  Reynolds  number. 

The  constant,  n,  can  be  calculated  by  the  following  equations 
given  by  Richa-dson  and  Zaki  (4). 

n  =  4.65  +   I9.5  <3/D 

n  -  (4.34  -h  17.5  d/D)  Ngg.03 

n  -=  (4.45  +18  d/D)  N"^'1 

n  *  4.45  »£♦* 

n  a  2.39 


NRe<  0.2 

(28) 

0.2<NRe<l 

(29) 

1    <    N_     <200 
Re 

(30) 

200    <    NRe<500 

(3D 

5°°   <    NRe 
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In  Table  3  all  of  the  experimentally  determined  time  constants 
are  presented  along  with  those  calculated  by  Equation  (1?).   The 
percentage  deviation  of  theoretical  and  calculated  tine  constants 
was  calculated  by 


$A   =  (Texp  -  Tcal)100 

T 

exp 


(33) 


In  Table  3  experiments  denoted  by  SU  were  for  a  step-up  in  fluid- 
izing-  velocity,  and  those  denoted  by  SD,  were  for  a  step-down. 
In  order  to  obtain  some  idea  of  tbe  applicability  of  the 
linearized  model,  it  was  necessary  to  devise  an  arbitrary  test. 
A  time  interval,  t-  nnt  was  defined  in  such  a  manner  that  the 
expression 


(h~hss ) 


-  (h~hss) 


exp 


cal 


< 


0.02 


(3*0 


was  valid  within  the  interval.  The  subscript,  exp  and  cal  in  ex- 
pression (3*0  refer  to  experimental  values  and  values  obtained  by 
the  linearized  model  respectively.  The  values  of  tQ  Q2  are  tabu- 
lated for  each  experimental  run  in  Table  ^. 
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Table  3.      St  en 

input  analysis* 

Exp.    No. 

nssv 

.)   :    e 

ss 

• 

A(cm. ) 

:    y(s*0.  ) 

:    i«(sec. 

)    !        %A 

• 

• 
• 

• 
• 

ex   . 

• 

.          Cal. 

• 
• 

SET   I   - 

70-80  mesh 

particles 

,    Bed  wt. 

400  gms, 

'       8       2,5 

g./cc. 

SU-2 

14.14 

0.658 

10.7 

22.2 

18.2 

18.0 

SU-3 

14.14 

0.658 

11.8 

21.6 

18.2 

15.7 

SU-4 

It,  14 

0.658 

2.3 

20.0 

18.2 

10.1 

SD-5 

33.44 

0.763 

19.2 

31.5 

38.1 

-20.9 

SD-6 

33.44 

0.763 

19.2 

31.8 

38.1 

-19.8 

SET   II 

-  60-70  mesh  particl 

M|    Bed  wt.        400   gms.      s 

2.5  g./cc. 

su-7 

16.2 

0.513 

33.5 

31.6 

30.1 

4.7 

SU-8 

22.8 

0.654 

6.6 

20.5 

21.0 

-2.4 

SU-9 

23.6 

0.664 

8.4 

20.9 

23.5 

-12.5 

SU-10 

16.1 

0.509 

7.9 

23.5 

20.6 

12.5 

SU-11 

15.9 

0.503 

31.5 

32.0 

32.3 

-0.9 

SD-12 

47.1 

0.832 

30.9 

40.0 

47.4 

-13.5 

SU-13 

14.8 

0.466 

9.2 

23.0 

23.1 

-0.4 

SD-14 

24.3 

0.675 

9.0 

20.0 

22.5 

-12.5 

SU-15 

15.6 

O.493 

18.7 

27.0 

27.2 

-0.7 

SD-16 

34.7 

0.772 

18.7 

28.5 

31.4 

-10.1 

SU-17 

15.7 

0.^97 

14.4 

tJuo 

24.6 

-2.5 

SD-18 

30.1 

0.738 

14.0 

24.0 

26.8 

-11.7 

SU-19 

16.7 

0.527 

8.4 

20.0 

19.6 

2.0 

SD-20 

25.1 

0.685 

7.5 

20.7 

21.8 

-5.3 

SU-21 

17.1 

0.538 

3.7 

22.0 

19.4 

7.3 

SB-?? 

20.8 

0.620 

3.7 

18.5 

19.6 

-5.9 

SU-23 

20.9 

0.621 

29.0 

37.0 

25.5 

31.1 

SD-24 

50.5 

0.843 

30.2 

41.5 

51.5 

-24. 1 

SU-2  5 

20.5 

0.615 

8*2 

19.0 

20.1 

-5.3 

SD-26 

28.7 

0*7*3 

8.2 

24.0 

25.1 

-4.6 

SU-2  7 

18.7 

0.577 

16.6 

°5.0 

24*2 

3.2 

SD-28 

35.4 

0.777 

15.4 

36.2 

35.1 

3.4 

SD-30 

^6.0 

0.859 

34.1 

56.0 

60.5 

-7.1 

SU-31 

25.6 

O.69I 

6.6 

20.5 

22.7 

•10*7 

SD-32 

32.? 

0.755 

6.7 

27.0 

28*4 

-5.2 

su-3  3 

23.1 

O.658 

6.1 

18.7 

20.6 

-10.2 

SD-34 

29.2 

0.730 

6.4 

22.5 

24.9 

-10.6 

su-35 

21.5 

0.632 

5.9 

19.8 

18.0 

9.1 

SD-36 

27.4 

0.712 

6.1 

21.0 

23.5 

-11.9 

su-37 

21.2 

0.627 

'.2 

of  Input 

18.5 
velocity; 

19.7 

SD  sign! 

-6.5 
fles  a 

*su  ■lgnlfie 

step-down. 

1   sto-p-up 
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Table 

3.      (contj 

Exp.    No. 

hss (cm.  ) 

• 

£"SS 

:    A(cm.) 

'      T(rec.)    : 

T(sec. ) 

:     %& 

• 

• 

• 

exp. 

cal. 

• 

• 

• 

I 

•                                                   • 

• 

SU-38 

20.0 

0.605 

9.6 

20.5 

19.4 

5.4 

SU-39 

17.3 

0.543 

9.9 

21.0 

20.3 

3.3 

SU-40 

22.8 

0.654 

6.5 

24.5 

21.3 

13.1 

SU-42 

23.* 

0.662 

12.3 

27.5 

30.5 

-9.1 

SU-^3 

24.1 

0.668 

11.9 

25.5 

24.2 

5.1 

SU-44 

18.1 

0.563 

9.6 

24.5 

?5.6 

-4.5 

SU-45 

17.5 

0.549 

6.5 

22.5 

24.4 

-8.5 

SU-46 

17.5 

0.549 

6.5 

23.5 

24.4 

-3.8 

SU-47 

29.1 

0.729 

9.6 

31.1 

28.8 

7.4 

SU-48 

29.1 

0.720 

9.4 

28.0 

28.8 

-2.9 

SET   III   - 

•  40-45  mei 

sh  parti 

cles,   Bed 

wt.      400  gmt 

5.             g       2,5 

g./cc. 

SU-49 

19.1 

0.587 

4.5 

8.5 

10.2 

-20.0 

SD-50 

23.6 

0.666 

4.0 

13.6 

12.3 

9.5 

SU-51 

19.1 

0.586 

5.1 

8.6 

9.1 

-17.4 

SD-52 

24.2 

0.673 

4.9 

11.8 

U.7 

0.8 

SU-53 

18.4 

0.571 

17.0 

11.0 

10.5 

4.5 

SD-54 

35.4 

0.777 

16.9 

19.4 

18.5 

4.6 

SU-55 

19.2 

0.586 

20.6 

14.0 

10.7 

23.5 

SD-56 

39.8 

0.801 

20.5 

19.5 

21.8 

-6.7 

SU-57 

26.2 

0.698 

3.9 

11.5 

13.0 

-13.1 

SD-58 

30.1 

0.738 

3.6 

11.5 

15.1 

-31.3 

SU-59 

?5.3 

0.687 

5.4 

10.8 

12.6 

-16.6 

SD-60 

30.6 

0.742 

5.7 

12.9 

15.2 

-17.8 

SU-61 

23.8 

O.667 

8.2 

10.5 

12.1 

-15.3 

SD-62 

32.0 

0.752 

8.3 

13.9 

16.1 

-15.8 

SU-63 

21.3 

0.628 

7.2 

10.3 

10.8 

-4.8 

SD-64 

28.4 

0.722 

7.1 

13.3 

13.7 

-3.0 

SU-65 

21.1 

0.625 

0.584 

14.3 

11.2 

10.7 

4.5 

SU-67 

19.0 

18.0 

12.0 

10.6 

11.6 

SD-68 

36.1 

0.780 

17.6 

15.0 

18.3 

-22.0 

SU-69 

I6.5 

0.520 

7.2 

6.5 

7.3 

-12.3 

SU-71 

15.6 

0.493 

10.1 

10.0 

10.8 

-8.0 

SU-73 

26.0 

0.685 

6.7 

13.5 

13.2 

2.2 

SD-74 

32.7 

0.758 

6.3 

I8.3 

16.4 

10.4 

SU-75 

16.1 

0.509 

23.8 

12.4 

10.6 

14.5 

SU-77 

32.8 

0.793 

16.1 

19.5 

20.1 

-3.1 

SU-79 

36.1 

0.781 

4.7 

18.4 

18.9 

-2.7 
-21.2 

SD-80 

40.8 

0.806 

3.8 

18.4 

22.3 

3« 

Table 

3,   (cont.) 

Exp.  No. 

:  hss(cm«) 

e  ss 

A(cm.)  : 

ip(sec. )  t 

>r(sec. ) 

:  %a 

• 
• 

•  • 

•  • 

• 
• 

exp.   . 

cal. 

• 
1 

SU-81 

26.2 

0.698 

4.8 

17.5 

13.0 

25.7 

SU-82 

18.3 

O.567 

7.7 

10.2 

9.8 

3.9 

SU-83 

18.5 

0.572 

6.5 

10.3 

10.0 

2.9 

SU-84 

22.6 

0.650 

11.6 

12.5 

11.7 

6.4 

SU-85 

24.3 

0.674 

11.6 

11.4 

12.1 

-6.1 

SU-86 

24.0 

O.670 

3.4 

11.5 

11.8 

-2.6 

SU-88 

20.6 

0.616 

18.4 

15.0 

12.9 

14.0 

SU-89 

20.7 

0.618 

21.3 

14.5 

13.1 

9.6 

SU-90 

21.0 

0.623 

15.3 

11.4 

11.6 

-1.8 

SU-91 

21.2 

0.627 

16.0 

11.4 

11.5 

-0.9 

SU-92 

21.4 

0.630 

4.9 

11.0 

11.6 

-5.5 

SU-94 

21.3 

0.628 

10.7 

13.9 

11.6 

16.5 

SU-95 

22.7 

0.652 

11.8 

14.5 

12.0 

17.2 

SU-96 

22.6 

0.650 

14.6 

15.6 

12.3 

14.7 

SU-97 

23.3 

0.660 

17.1 

14.8 

13.3 

10.1 

SU-98 

23.3 

0.660 

21.0 

12.0 

13.3 

-10.8 

SU-99 

23.9 

0.669 

21.6 

12.9 

12.0 

7.5 

SU-100 

24.0 

0.670 

15.0 

16.0 

12.3 

23.0 

SU-101 

24.4 

0.675 

11.4 

15.0 

12.8 

14.6 

SU-102 

24,2 

0.673 

16.4 

15.2 

12.3 

19.1 

SET  IV  - 

25-30  mesh 

particles 

,  Bed  wt. 

200  gms. 

»   s   ** 

08  g./cc. 

SU-103 

27.2 

0.660 

25.4 

60.0 

67.1 

-11.8 

SU-104 

23.1 

0.599 

24.1 

41.0 

48.7 

-18.8 

SU-105 

20.4 

O.547 

26.1 

57.0 

59.9 

-5.1 

SU-106 

20.2 

0.543 

11.6 

46.5 

53.2 

-14.4 

SU-107 

22.5 

0.590 

12.3 

46.5 

51.2 

-10.1 

SU-108 

23.8 

0.612 

18.8 

49.5 

51.9 

-4.8 

SU-109 

19.9 

0.535 

31.1 

44.0 

48.4 

-10.0 

SU-110 

19.6 

O.528 

16.0 

49.0 

48.2 

1.6 

SU-111 

19.8 

0.533 

8.2 

52.4 

48.4 

7.6 

SD-112 

26.6 

0.652 

6.6 

63.7 

58.6 

8.0 

SU-113 

20.1 

0.539 

12.7 

51.5 

48.6 

5.6 

SD-114 

32.8 

0.718 

12.5 

68.7 

71.5 

-4.7 

SU-115 

23.3 

O.603 

17.7 

50.5 

51.2 

-1.4 

SD-116 

41.3 

0.776 

17.3 

71.0 

78.8 

-11.0 

SU-117 

25.3 

0.635 

10.9 

44.0 

55.9 

-27.1 

SD-118 

36.2 

0.744 

7*2 

56.3 

74.6 

-32.5 

• 
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Table  k, 

Values 

of  t0#0?  in 

seconds. 

Exp.    No.* 

t0.02 

Exr>.    No 

•          fc0.02 

Exp.    No. 

t0.02 

SU-2 

11 

88*28 

6 

SD-68 

2 

SU-3 

SD-30 

9 

SU-71 

10 

SU-4 

55 

80-31 

27 

SU-73 

4 

SD-5 

12 

SD-34 

,1 

SD-7^ 

6 

SD-6 

7 

SC-35 

SU-75 

28 

SU-7 

75 

SD-36 

7 

SU-77 

30 

SU-8 

17 

SU-49 

4 

3U-79 

6 

SU-9 

18 

su-50 

6 

SD-80 

3 

SU-11 

88 

KJ-51 

20 

SU-82 

6 

SD-12 

8 

SD-52 

7 

SU-8  5 

23 

su-13 

I 

su-53 

36 

SU-89 

19 

SD-14 

SD-54 

6 

SU-106 

35 

SU-15 

80 

su-55 

37 

SU-107 

33 

SD-16 

7 

SD-56 

2 

SD-112 

13 

SU-17 

su-57 

4 

SU-113 

26 

SD-18 

8 

SD-58 

3 

SD-114 

18 

SU-19 

35 

SU-59 

3 

SU-115 

85 

SD-20 

P 

SD-60 

4 

SD-116 

16 

SU-21 

SU-61 

5 

SU-117 

70 

SD-22 

12 

SD-62 

3 

SDull8 

11 

SU-23 

150 

SU-63 

12 

SD-24 

a 

SD-64 

2 

SU-25 

SU-65 

17 

SD-26 

7 

SU-67 

27 

*SU  signifies  a 

step-up  of 

input  velocity;    SD  signifi 

es  a 

step-down. 
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Response  to  Pulse  Input 

A  considerable  amount  of  oulse  data  was  taken  In  the  manner 
described  In  the  experimental  section.   Figure  10  shows  the  re- 
sponse observed  for  a  representative  experimental  run.   A  dis- 
cussion of  the  deviation  of  the  experimental  response  from  the 
model  can  be  found  in  the  section,  "Deviations  from  the  Model." 
The  purpose  for  obtaining  data  in  this  form  was  to  convert  it  to 
frequency  response  data  by  an  approximate  Fourier  transform  tech- 
nique.  Hougen  and  Walsh  (6)  have  discussed  "pulse  testing"  as  a 
convenient  method  of  obtaining  frequency  response  data.   A  pulse 
is  described  as  a  function  which  differs  from  zero  for  a  finite 
time  only.   The  form  of  the  pulse  applicable  to  the  pulse  testing- 
method  is  arbitrary  and  in  these  experiments  consisted  of  a  com- 
bination of  two  step  changes  as  previously  described. 

A  time  function,  f(t),  is  composed  of  an  infinite  number  of 
harmonic  components.   The  resolution  of  these  components  is  ac- 
complished by  evaluation  of  the  Fourier  transform 


F(i<o  )  -  f  f(t)  e~lu;t 


dt  (35) 


A  method  and  tables  for  evaluation  of  this  integral  are  given 
by  Huss  and  Donepan  (?).   The  method  presented  consisted  of  repre- 
senting f(t)  by  a  staircase  function  consisting  of  n  equal  inter- 
vals.  This  procedure  is  shown  in  Figure  11.   The  Fourier  trans- 
form in  terms  of  the  transfer  'unction  becomes 
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F(i^>  )  •»       H   F„(l  oo  ) 


(36) 


where 


P  (i^  )  = 


^.i 


rn  e-1  "  *  dt 


The  re^  part  of  Equation  (37)  Is 


B 


n 


(l^o  )]=      r     2   Bln(Ai.6A)  )   cos        (2n-l)  4JL 
J  n  co  p  o 


to 


(37) 


(38) 


and  the  Imaginary  part   is 
In  [P(i^  ) 


-     r^  ~  sin(A£.ou)   sin 


(2n-l)  A±oo 
2 


(39) 


where  r     is   the  height  of  the   staircase  function  and   4  t   is   the 
time  interval. 

The   real   part  of   the  total   Fourier  transform  of  the  function 
is  the   sum  of   the   real  parts  of  the   individual   intervals. 


R 


F(ico   >    =.      H    \    [p(1uj   )J 


(40) 


and  the   imaginary  part   is  the  sura  of  the   imaginary  parts  of  the 
intervals 


Letting-  z-   -~tco    in  Equation    (38),    Equation    (40)   becomes 


(41) 


[f(1^  jj~    ^it£ra 


R        F(i^i 


sin   |   cos(2n-l)z 


(42) 
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and   from    Equation    (39)   and   Equation    (41) 


I     [p(i*o    )]     «     -    <4t   £  rTi?ln   7   sln(2n-Dz  (43) 


The  modulus  of  the  Fourier  transform   is 


A(^   )  =     J  %JL   [p(icO  )g     2_^     £    [F(i^>  j/j     2  (44) 

and  the  phase  angle 

**<*)  ==     -   tan-1   *  F(i^   ?  (45) 

HP(ioo) 

The  evaluation  of  Equations  (42)  and  (43)  is  simplified  by  use  of 

tables  in  which  values  of  the  functions  sin  z  cos(2n-l)z  and 

z 

sin  z  sln(?n-l)z  are  tabulated  for  values  of  n  and  z   (?).   The 
z 

amplitude  ratio  can  be  determined  as 

A.R.  =r        A(C°    *    Mtittttl  (46) 

A(<°  }    in^ut 

and  the  nhase  lag  is 

^  ^inr,ut  -   ^output  <47> 

This  calculation  was  carried  out  for  experimental  pulse  in- 
puts of  from  15  to  30  second  duration.   One  of  the  better  response 
curves  obtained  is  presented  in  Figure  12. 

The  pulse  testing  technique  requires  that  the  system  be 
linear.   Application  to  a  non-linear  system  will  not  give  meaning- 
ful results  unless  the  disturbance  considered  is  small  enough  for 
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the  system  to  behave  line-rly.   When  the  pulse  data  were  being 
taken  it  was  not  known  wh3ther  this  condition  was  being  fulfilled 
or  not.   However,  due  to  the  rather  odd  shape  of  the  response 
curve,  particularly  that  for  phase  lag,  it  MM  considered  necessary 
to  measure  the  actual  resoonse  to  a  sinusoidal  input. 

Response  to  Sinusoidal  Inr>ut 

Two  types  of  data  were  taken  on  the  response  to  a  sinusoidal 
input.  The  first  was  bed  height  versus  time  data.  This  was  done 
in  order  to  determine  the  degree  of  distortion  in  the  output.  A 
linear  system  will  produce  a  response  to  a  sinusoidal  input  which 
iB  also  sinusoidal.  Non-linearities  will,  however,  tend  to  cause 
deviations  from  the  sine  w^ve  form.  No  attempt  was  made  to  quan- 
titatively analyze  the  distortion,  but  a  graph  illustrating  this 
behavior  is  ^resented  in  Figure  13. 

The  second  type  of  data  taken  was  used  in  the  determination 
of  amplitude  ratio  and  phase  lag  relationships.  Amplitude  ratio 
Is  defined  by 

A.R.  -   amplitude  cf  output  (48) 

amplitude  of  input 

The  output  amplitude  was  determined  as  the  difference  between 
the  maximum  and  minimum  in  bed  height  expansion  divided  by  two. 
The  input  amplitude  was  a  sinusoidally  varying  flow  rate.   For 
use  in  computation,  this  had  to  be  converted  into  the  same  units 
^s  th^s  outrut  amplitude.   This  was  accomplished  by  UBe  of  the  stead- 
y-st^te  curve  relating  bed  height  and  flow  rate  for  the  system 
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being  studied.   The  procedure  used  is  illustrated  in  Figure  14. 
In  this  figure,  the  input  flow  amplitude  is  0.32  cm. /sec.  and 
the  steady-state  flow,  2.16  cm. /sec.   Thus,  the  maximum  and  mini- 
mum in  flow  are  2.48  cm. /sec.  and  1.84  cm. /sec.  respectively. 
The  steady-state  heights  corresponding  to  these  flows  are  37*0 
cm.  and  25.0  cm.   This  gives  an  amplitude  in  terms  of  bed  height 

expansion  of  37tO-?'?>0   6.0  cm.   The  mechanism  by  which  the  me- 

2 

chanical   sine  wave  generator  r>roduces  a   given  amplitude   in  flow 
can  be  seen  by  referring  to  Figure  5«      From  the   ex-ores  si  on 

V    »     Acy  VL  sln/f  (49) 

it   is   obvious   that   the  flow,    V,    is   a   maximum  when    /&  -=•      ^2  and 
a  minimum  when   /?  —      3  tt  /2.      The  values  of  these  extremes  are 

V#  »    ±     Acy  VL  <50) 

The   tangential  velocity,   V^,    is   given  by 

VL    =      LuJ  (51) 

where  L  is  the  length  of  the  lev  r  arm  and  <^0    is  the  frequency. 
Therefore,  from  Equation  (50) 

V*=  *   Acy  L^  '  (52) 

or  in  terras  of  superficial  velocity 

ool  Acol 

where  <j>  *   is    J.he  superficial   liquid  velocity  amplitude  and  Acol   is 
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the  cross-sectional  area   of   the  column.      Thus  various   input  ampli- 
tudes were  obtained  by  varying  the  values   of   L  and  <-^  . 

The  experimental  phase  lag  was  determined  by  obtaining  the 
time  interval  between  the  reversal  in  piston  motion  and  the  re- 
versal in  bed  height  motion.  The  phase  lag  was  then  determined 
as 

tan   (90°  -<)   «      k  ZUAMl       0°<  «<  90°  (54) 

2  T 

where  °<    is  the  T)hase  lag  and  A  t  is  the  tima  interval  between  re- 
versal in  the  direction  of  nlston  motion  end   the  reversal  of  bed 
height  motion.   This  method  was  very  limited  since   t  is  quite 
small  at  high  frequencies  and  could  not  be  measured  by  observation. 
Another  limitation  lies  in  the  fact  that  the  wave  form  of  the  out- 
put is  distorted  as  illustrated  in  Figure  13* 

Amplitude  and  phase  lag  data  are  shown  plotted  against  fre- 
quency in  Figure  15  for  the  6  to  7  mesh  glass  particles.  In  this 
series  of  experiments  no  attempt  was  made  to  hold  the  input  ampli- 
tude constant  for  each  experiment.  The  input  amplitude  range  is 
from  4.05  cm. /sec.  to  0.494  cm. /sec.  or  from  10.4  cm.  to  1.4  cm. 
in  terms  of  bed  height  expansion.  These  data  are  shown  compared 
with  the  response  calculated  by  the  linearized  model* 

DISCUSSION  OF  RESULTS 

At  this  point  it  would  be  desirable  to  emphasize  the  nurpose 
of  the  work.   The  basic  desire  was  to  develon  a  linearized  model 
to  describe  the  dynamic  behavior  of  liquid-solid  fluidized  beds 
and  to  show  by  experiment  that  the  model  is  sufficiently  accurate 


1.0 
0.8 

— ghtta      l      1      •      '      •      •      •      '      ' 

No 

o 

o  0.6 

^«^S 

0) 

_                                                                                                                  ^^"vv^^     — 

.?  0.4 

Glass    bead      6-7  mesh 

■ 

a 

E 
< 

O    Experimental    data 
—  Linearized    model 

0.2 
0 

T  s  1.5  sec 

_  10 

0) 

£ 

«30 

o 

^o^o 

""  50 

—                                                                                                                                            ^"^-~v^^                                                                               " 

0) 
(0 

o 
1  70 

—                                                                                                                                      ~"— «-^,,~" 

90, 

1      1       1      1      1       1      1       1      1       1       1       1 

D          0.2         0.4         0.6         0.8          1.0           1.2 

(£      (rad/sec) 

Fic 

).  15.     Frequency     response    of     liquid- 

solid     fluidized     bed     for    6-7 

mesh    glass    beads. 

1 

50 


for  control  purposes.   Many  non-linear  flow  processes  are  treated 
in  this  manner  for  the  purpose  cf  ar-ly-'-  **4   synthesizing  pro- 
cess  control  systems.   Any  fluid  flow  or  heat  transfer  process 
will  exhibit  non-linearities  due  to  the  dependence  of  transport 
properties  on  the  values  of  such  intensive  variables  as  tempera- 
ture and  pressure.   For  example  in  hent  transfer  the  conductivity 
of  a  material  is  dependent  on  the  temperature,  and  in  fluid  flow 
the  resistance  to  flow  is  usually  dependent  on  Reynolds  number, 
which  is  a  function  of  fluid  density,  velocity,  and  viscosity. 
In  some  instances  it  is  possible  to  very  nearly  approximate  system 
behavior  by  a  linear  differential  equation;  while  in  other  cases 
this  might  not  be  true.   The  accuracy  with  which  a  system  may  be 
approximated  is  usually  dependent  on  the  relationship  between  the 
steady-state  behavior  of  the  system  and  the  independent  variables. 
If  this  relationship  were  linear,  we  could  employ  a  linear  differ- 
ential equation  to  describe  the  dynamic  behavior  of  the  system 
exactly.   Therefore,  we  should  examine  the  steady-state  relation- 
ship between  fluidizing  velocity  and  porosity  or  bed  height. 
This  relationship  is  given  by  Richardson,  et.  al.  (k)   as 

Since  norosity  is  the  independent  variable  in  this  case  it  would 
be  well  to  rewrite  this  relationship  as 


6»«  k<#/n 


where 
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*  -  !/Us 


1/n 


The  terminal  falling  velocity,  Us,  will  remain  constant  when  the 
fluid  viscosity  and  density  are  constant.   Therefore,  non-linearity 
of  this  curve  depends  on  the  value  of  1/n.   Prom  Equations  (28)  to 
(3D  it  is  seen  that  in  the  usual  range  of  operation,  the  value  of 
n  ranges  from  2.39  to  4.65.   Values  of  n  within  this  range  do  not 
depart  sufficiently  from  one  to  make  linearization  impractical. 
There  is,  however,  another  factor  which  makes  fluidizatlon 
very  difficult  to  analyze  from  a  dynamic  standpoint.   This  involves 
the  mechanism  of  bed  expansion  and  contraction.   As  shown  by  Slis, 
et.  al.,  (5),  the  mechanism  for  expansion  is  different  from  that 
for  contraction.   Because  of  this  difference  between  expansion 
and  contraction,  the  response  to  a  decreasing  flow  rate  differs 
fror^  that  to  an  increasing  flow  rate.   Therefore  it  is  Impossible 
to  describe  the  dynamic  behevior  of  a  fluldized  bed  accurately  by 
one  rodel.   This  is  true  whether  the  model  is  linear  or  non-linear. 
It  should,  however,  be  possible  to  describe  the  dynamic  behavior 
for  small  fluctuations  by  a  single,  linear,  model  and  this  is  what 
has  been  attempted.   The  necessity  for  this  type  of  representa- 
tion is  obvious  when  it  is  considered  that  the  deviations  which 
might  occur  in  an  actual  process  are  random. 

Effect  of  Changes  in  Temperature  on  the  Time  Constant,  T 

Since  many  processes  may  be  subject  to  changes  in  fluid  tem- 
perature it  is  necessary  to  consider  such  effects  on  the  time  con- 
stant of  the  r^roces?.   The  expression  for  the  time  constant,  T, 
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is  given  as 

h 


T  - 


ss 


r,Us  G»-»  (1-y 


or  the  equivalent  expression 


T  =     ss  ^  ss 


At  a  constant  bed  height  and  porosity  the  dependence  of  T  on  tem- 
perature is  seen  to  involve  the  changes  of  n  and  <£  or  n  aild  Us. 
Since  n  is  an  empirical  function  of  Reynolds  number  and  Ug  is  a 
function  of  liquid  density  and  viscosity  it  would  be  useless  to 
try  to  develcoe  an  analytical  expression  for  this  effect.   In 
Table  5,  however,  values  of  T  are  given  as  a  function  of  tempera- 
ture for  a  system  consisting  of  water  and  40  to  ^5  mesh  glass 
beads.   Table  5  shows  that  a  considerable  change  in  time  constant, 
T,  might  be  expected  when  temperature  variations  of  ±.  50°F  or 
greater  are  encountered.   For  smaller  changes  in  temperature  this 
effect  would  probably  not  be  important. 

Deviations  from  the  Model 

For  the  response  to  a  stet>  lnr>ut,  there  are  two  tynes  of  de- 
viations to  be  considered.   The  first  involves  the  deviations  of 
the  time  constants  calculated  by  the  linearized  model  from  the 
experimental  time  constants.   These  deviations,  expressed  in  per- 
centage, have  been  tabulated  for  e^ch  experimental  run  in  Table  3. 
By  considering  both  positive  and  negative  deviations,  the  average 
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Table   5»      Dependence  of  T  on  Temperature. 

Particles:      40  -  45  mesh  glass  beads    (  />  2.5  g./cc.) 

Fluid:      Water,   hQJ-    25  cm.     £  oa  -  0.600. 

TeroTD.    (op)                   Us    (cm. /sec.)                     n  T    (sec.) 

50                                       4.00                            3.46  15.86 

75                                    4.4-5                         3.33  13.83 

100                                    4.98                         3.21  12.05 

150                                    6.24                         3.00  9.28 

200                                       7.53                            2.82  7.48 


5^ 


deviation  for  ill  experimental  runs  was  -2.1  %.      Since  no  systema- 
tic variation  in  these  deviations  could  be  related  to  any  of  the 
operational  variables  such  as  bed  height,  porosity,  or  size  of 
step,  it  must  be  assumed  that  they  are  the  result  of  inaccuracies 
in  calculation  and  the  inherent  inaccuracy  in  the  graphical  method 
of  obtaining  experimental  time  constants.   On  the  other  hand,  the 
time  interval  over  which  the  experimental  response  follows  the 
linearized  model  is  definitely  a  function  of  the  operating  condi- 
tions.  This  has  been  shown  in  a  previous  section  where  t0>02  was 
tabulated  for  each  run.   It  is  evident  that  in  general  the  re- 
sponse to  a  step  up  is  described  accurately  over  a  much  greater 
time  interval  than  the  response  to  a  step  down.   Because  of  this 
it  is  necessary  to  consider  the  resoonse  to  a  sten  down  in  flow 
in  order  to  determine  the  limitations  on  the  applicability  of  the 
linear  model. 

Slis,  et.  al.  (5)  have  shown  that  the  response  to  a  step  down 
in  f  lev.  can  be  considered  a  special  case  of  sedimentation.   There- 
fore, we  expect  a  constant  bed  height  velocity  resulting  in  a 
straight  line  response  of  bed  height  to  a  step-down  in  fluldlzlng 
velocity.   This  response  can  be  represented  by 

*  »  ;<4>X*  *ss)  t  +  hss  (55) 

where  $  t  ll  the  superficial  llauid  velocity  after  the  steo  innut. 
The  linearized  model,  however,  gives  a  response  of 

h  -  A  (l-e-^)  +  hBa  (56) 


55 

rearranging   (55)   and    (56)   and  taking  their 

difference  gives 

(5±tt)            -    (!2±£&)        -     kt  - 

tW**J 

(57) 

hyp               *       cal 

where  k   -i^l"     «b).      Here  the   subscript 

(hyo)    stands  for  a 

hy- 

oothetical  response  based  on  Equation   (55) 

and    (cal)    for  the 

re- 

spons 

*e  calculated   from  the  linearized  model.      Comparing   Equa 

tlon 

C57) 

with   y.quatl~>n    (3M    it   is   seen  that 

t_).0? 

k  tlt0?  -  (l-e      ■       ;  -  0.02 

(58) 

This 

can  be  rewritten  as 

t0.0? 
e                T        -«-        k  t%02    « 

1.02 

(59) 

UtlM   (59)   can  be  further  simplified  by  noting  that 

A^>~ 

.  «p  from  Equation    (13).      Since  k      A-jX-  , 

we  can  write 

k  -    hH~  i- 

T 

(50) 

■ai 

(59)   becomes 

*-- H1    -  Sbii  -  1.02 

1 

(£1) 

Now 

Eoustlon   (61)   is  a   single  valued  function  of   (to  q?)    so 

that 

there  will   be  only  one  value  of    {^'      )  which  will   satisfy  it. 

This 

is  graphically  found  to  be   (tQ  Q2/"r^ 

—      0.205.      '^be  actual 
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values  of  tor02   from  Tables  3  and  k   vary  from  0.1  to  0.35  with  a 

T 
mean  of  0*22,   Part  of  this  error  Is  undoubtedly  due  to  the  trial 

and  error  method  of  obtaining  the  t0#02  and  the  remainder  to  the 
inaccuracy  of  Equation  (61).   Equation  (61)  would,  however,  be  an 
acceptable  means  of  estimating  the  deviation  to  be  expected  when 
applying  the  linearized  model  to  a  step-down  in  fluidizing  veloc- 
ity.  It  should  be  emphasized  that  the  value  0.02  was  chosen 
arbitrarily  and  that  any  desired  value  may  be  picked.   This  is 
shown  by  writing  Equation  (61)  in  the  more  general  form 

e-t/T+   |  *f  (t/T)  (62) 

T 

where 

h— h  h— h„„  , ,    . 

f<t/T)^  i+ <— fit)        -    (— 7^)  <63) 

A   hyp        A   cal 

Therefore  any  desired  value  of   (  "  ss)    -  (  "  ss)    can  De 

A  hyp       A   cal 
chosen  and  the  corresponding  value  of  (t/T)  calculated.   To  facil- 
itate this  determination,  a  plot  of  f (t/T)  versus  t/T  is  given  in 
Figure  16. 

The  response  to  a  step-up  in  fluidizing  velocity  Is  generally 
approximated  very  closely  by  the  linear  model  over  a  long  time  in- 
terval.  An  Interesting  phenomenon  which  was  observed  concerning 
the  response  to  a  step-up  was  the  development  of  an  irregularity 
in  the  response  curve  when  the  input  was  quite  large.   It  is  not 
known  whether  this  was  due  to  some  experimental  error  or  whether 
it  is  inherent  in  the  response  to  a  step-up  in  fluidizing  velocity. 
The  development  of  this  irregularity  is  shown  in  Figures  17a,  17b 
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and  l?c.   As  can  be  seen,  the  linearized  model  will  adequately 
represent  the  response  over  tMs  irregularity,  but  deviates  there- 
after.  Although  this  behavior  is  interesting,  it  probably  would 
not  ordinarily  be  of  practical  interest  since  the  irregularity 
was  never  observed  to  develop  except  when  the  disturbance  was  of 
a  magnitude  equivalent  to  the  steady-st^te  height.   In  process 
control  applications  it  is  improbable  th^t  disturbances  of  this 
magnitude  would  be  encountered. 

The  same  non-linearities  which  caused  deviations  in  the  re- 
sponse to  a  step  input  were  -^lso  found  significant  in  the  fre- 
quency response  analysis. 

Insufficient  frequency  response  dat^  has  been  obtained  to 
permit  a  complete  analysis  of  the  deviations  of  this  response  from 
the  model  and  about  all  that  is  possible  is  a  qualitative  des- 
cription.  Figure  13  shewed  that  a  small  amount  of  distortion  in 
the  output  wave  form  was  obtained.   This  is  reasonable  and  can 
be  attributed  chiefly  to  two  sources.   First,  the  input  in  terms 
of  flow  rate  is  a  true  sine  wave  (subject  to  experimental  error), 
but  the  input  in  terms  of  equivalent  bed  height  is  not.   This  is 
because  of  the  non-linear  nature  of  the  steady-stnte  relation  be- 
tween bed  height  and  fluidizing  velocity  which  has  been  clearly 
illustrated  in  Figure  14.   Therefore,  since  the  input  is  not  ex- 
actly sinusoidal,  we  could  not  expect  the  output  to  be.   This 
source  of  distortion  is,  of  course,  dependent  on  the  input  ampli- 
tude and  becomes  negligible  as  the  input  amplitude  approaches 
zero.   The  second  source  of  output  distortion  is  the  difference 
in  the  mechanisms  for  expansion  and  contraction  of  the  bed.   This 
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effect  could  be  quite  important  and  should  not  be  particularity 
related  to  the  input  amplitude.   Therefore  a  small  amount  of  dis- 
tortion can  be  expected  for  any  input  amplitude.   It  should  be 
noted  that  the  response  curve  presented  in  Figure  13  is  the  result 
of  rather  extreme  conditions.   The  input  amplitude  in  this  case 
is  14.4  cm.  or  approximately  0.5  h^.   Also,  the  manner  of  con- 
structing  Figure  13  tends  to  concentrate  all  of  the  distortion  in 
the  lower  half  of  the  curve.   The  reason  for  this  is  that  one  ex- 
perimental rioint  had  to  be  selected  as  a  reference  point  for  the 

remainder  of  the  curve.   This  was  chosen  as  the  point  at  (£^i£.) 

2tt 

equal  to  0.75.   In  the  actual  response  it  is  doubtful  that  all  of 
the  distortion  would  be  concentrated  in  one  area  of  the  curve  as 
it  is  shown.   The  author  feels,  however,  that  due  to  the  mechanism 
for  bed  contraction,  (the  bed  tends  to  contract  at  a  higher  rate 
than  it  expands  for  the  same  input),  more  distortion  would  be  ex- 
pected in  the  lower  half  of  the  cycle  than  in  the  upper  half  and 
Figure  13  was  constructed  in  such  a  manner  to  emphasize  this. 

Very  little  data  was  obtained  on  the  phase  lag  relationship. 
Experimental  limitations  of  the  procedure  as  described  and  the  dis- 
tortion in  the  output  wave  made  this  data  very  difficult  to  obtain 
and  interpret.   The  limited  amount  of  data  presented  was  obtained 
by  averaging  the  ohase  lag  at  the  top  of  the  cycle  with  that  at 
the  bottom.   This  procedure  generally  gave  results  in  fairly  good 
agreement  with  the  linearized  model.   One  noint  in  question  at 
this  time  concerns  the  limits  of  the  phase  lag  when  expressed  as 
a  function  of  frequency.   The  first  order  linear  model  predicts 
that  the  lag  should  be  zero  degrees  at  a  frequency  of  zero  and 
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approach  90  degrees  as  the  frequency  approaches  infinity.   This 
behavior  was  borne  out  over  the  entire  experimental  range  except 
at  the  maximum  frequency  obtainable  with  the  equipment.   At  this 
frequency,  approximately  0.92  rad./eec,  it  appeared  that  the 
phase  lag  was  Just  over  90  degrees.   This  behavior  did  not  mani- 
fest itself  to  a  sufficient  degree  M  this  frequency  that  any 
definite  conclusions  could  be  drawn  and  it  was  impossible  to  ob- 
tain higher  frequencys  with  the  equipment  as  constructed.   This 
would  be  an  interesting  area  for  investigation,  however,  since 
greater  lags  than  90°  would  suggest  the  propriety  of  models  of 
higher  order  than  one. 

The  amplitude  ratio-frequency  relationship  is  probably  the 
most  significant  tool  which  the  field  of  process  dynamics  has  to 
offer  for  use  in  the  design  and  evaluation  of  processes.   It  was 
found  that  the  first  order  linear  model  could  be  used  to  represent 
this  relationship  quite  accurately  for  reasonable  input  amplitudes 
{A   h  <.  0.5  hgg).   It  is  not  known  exactly  to  what  degree  of 
accuracy  the  model  could  be  expected  to  approximate  any  given 
system,  but  It  seems  likely  that  it  could  be  depended  on  to  yield 
results  satisfactory  for  preliminary  design  calculations.   In 
addition  to  the  results  given  in  Figure  15  for  6  to  7  mesh  glass 
beads  with  a  variable  input  amplitude,  additional  amplitude  ratio 
data  are  presented  in  Figure  18.   The  agreement  of  the  data  shown 
with  the  curve  calculated  using  Equation  (23)  is  very  good. 

A  final  point  on  the  deviations  encountered  in  the  response 
to  a  sinusoidal  input  is  concerned  with  the  "shift"  in  the  steady- 
state  height,  hss,  which  was  observed.   This  car  be  seen  in  Figure 
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13.   There,  the  midnoint  of  the  output  cycle  falls  about  0.7  cm. 
below  the  steady-state  height.   This  behavior  was  observed  in  all 
of  the  experiments  and  the  "shifts"  were  usually  downward.   From 
Figure  14  it  would  seem  that  any  change  in  hss  should  be  upward 
and  in  some  cases  a  slight  tendency  toward  this  was  observed.   In 
all  other  cases  it  seems  that  some  other  factors  outweigh  the 
effect  shown  in  Figure  14.   This  was  probably  caused  by  the  tend- 
ency for  the  bed  to  contract  more  rapidly  than  to  expand  as  men- 
tioned previously. 

Application  of  Pulse  Testing 
Method  to  Fluidizatlon 

As  has  been  Indicated,  an  attempt  was  made  to  obtain  frequency 
response  data  from  the  response  to  a  oulse  input.   An  example  of 
the  results  obtained  was  presented  in  Figure  11.   The  results  were 
generally  very  erratic  and  in  many  cases  the  amplitude  ratio  would 
not  approach  one  as  the  frequency  approached  zero.   Also  there 
was  a  tendency  for  the  analysis  to  yield  very  high  values  at  a 
frequency  around  0.4  rad./sec.   At  intermediate  frequencies, 
(0.1  -  0.3  rad./sec),  the  method  usually  gave  reasonable  results. 
The  poor  results  at  high  and  low  frequencies  could  be  caused  part- 
ly by  a  deficiency  in  the  harmonic  content  of  the  pulse.   The 
pulse  used  was  a  rectangular  one  and  the  response  to  this  pulse 
form  is  very  limited  in  harmonic  content.   It  would  be  advisable 
to  employ  a  more  Irregular  pulse  Input  in  any  further  work  along 
this  line.   The  non-linearity  of  the  system  also  contributes  very 
definitely  to  the  erratic  results.   It  must  be  concluded  that  pulse 
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testing  is  not  an  advisable  manner  in  which  to  obtain  frequency 
response  data  for  bed  height  change  in  fluidization  processes. 

Comparison  of  Results  with  Those 
of  Other  Investigations 

There  are  only  two  references  to  the  authors  knowledge  which 
present  any  results  on  the  study  of  transient  behavior  of  the 
fluldized  bed.   The  first  of  these  is  a  paper  by  Slis,  et.  al. 
(5)  in  which  the  response  to  a  step  input  in  fluldizing  velocity 
is  considered.   The  data  given  in  that  paper  is  not  of  such  a  na- 
ture that  a  direct  comparison  can  be  made  with  the  results  of  this 
present  work.   The  experimental  responses  observed  were  similar 
with  a  single  exception  which  will  be  discussed  under  the  topic 
of  "Effect  of  Particle  Size  Distributions." 

The  second  reference  is  a  technical  report  by  Martin  Nuclear 
Company  (8).   In  this  report  the  following  differential  equation 
was  derived  to  represent  the  transient  behavior  of  fluidized  bed 
height. 

d2L      g(/y^c)       6\   3^-2  £  ^^  ^  ^ 
dt2"     /%       e  €        At 

where  L  is  the  bed  height,  Z7  ^  the  solid  density,  yaQ   the  liquid 

density,  g  the  gravitational  constant,  uQ  the  superficial  liquid 

velocity,  €  the  porosity,  and  b  Is  a  constant  defined  by  the 

equation  6  BS  =•  du:?   . 

ss 
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No  experimental  data  was  presented,  but  Equation  (63)  MM 
programed  on  an  analog  computer  and  the  traces  for  various  inputs 
reproduced.   The  linear,  first  ordor  -nodel  was  used  to  calculate 
the  response  to  these  same  inputs  and  the  results  are  given  in 
Figure  19.   A  detailed  description  of  the  manner  in  which  the 
calculations  were  made  to  obtain  the  response  given  by  the 
linearized  model  is  presented  in  Appendix  B.   Figure  19  represents 
only  a  comparison  between  the  two  models  since  no  experimental 
data  is  available  with  which  to  compare  the  response  given  by 
Equation  (63).   The  curve  labeled  "Response  by  Equation  (63)" 
is  an  analog  response  which  was  given  in  reference  (8).   This  rep- 
resents an  interesting  application  of  the  linearized  model,  how- 
ever, since  the  operating  conditions  differ  widely  from  those  of 
the  experimental  conditions  existing  in  this  present  work.   This 
can  be  seen  by  referring  to  Appendix  B. 

Effect  of  Particle  Size  Distributions 

This  discussion  was  prompted  by  a  discrepancy  in  the  experi- 
mental results  of  this  work  and  the  results  presented  in  a  taper 
by  Slls,  et.  al.  (5).   They  reported  an  experimental  response  to 
a  step-up  in  flow  which  was  linear  for  the  first  part  of  the  re- 
sponse.  It  turned  out,  in  fact,  that  the  time  interval  over  which 
this  linear  response  was  observed  was  defined  by 

0  <  t  <  t* 

where 
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t»  •=-        e   sg  hsg  (64) 

neb      {US     ) 

'  ss     S8 

It  should  be  noted  th«t  this  equation  for  t*  Is  the  same  equation 
used  to  calculate  T.   In  this  present  experimental  work  the  time 
constant,  T,  was  usually  In  the  neighborhood  of  10-20  seconds.   It 
should  have  been  possible  to  detect  this  linear  relationship  quite 
easily,  but  It  was  not  observed.   The  response  to  a  step-down  was 
quite  definitely  linear  as  predicted  In  reference  (5).   This  a- 
greement  In  one  case  and  disagreement  In  another  related  case  is 
very  perplexing.   A  possible  explanation  Involves  the  effect  of 
particle  size  distributions  on  the  response  of  fluidized  beds  to 
upsets*   In  any  fluidization  n-ocess  there  will  be  a  certain  range 
of  particle  diameters.   If  all  of  the  particles  are  of  uniform 
density,  as  is  usually  the  case,  the  smaller  p-rtlcles  will  tend 
to  concentrate  -t   the  top  of  the  bed  and  the  larger  particles  at 
the  bottom.   Therefore  we  might  expect  a  somewhat  continuous  ax- 
ial distribution  of  particle  sizes.  For  particles  of  uniform 
density,  the  terminal  falling  velocity  will  increase  with  Increas- 
ing diameter.  This  will  also  cause  an  increase  in  Reynolds 
number  and  a  corresponding  decrease  in  the  empirical  constant,  n. 
However,  the  increase  in  U8  will  always  be  proportionally  greater 
than  th?  decrease  In  n.   deferring  to  Equation  (1?)  It  is  seen 
that  the  time  constant  will  decrease  with  increasing  particle 
diameter  due  to  the  effects  described  above  ^nd  larger  particles 
can  be  expected  to  show  a  faster  resnonse  to  upsets  than  smaller 
particles.   Now  if  the  bed  is  at  a  steady-state  condition  and  Is 
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subjected  to  a  step-ur  in  fluidizlng  velocity,  the  particles  at 
the  top  of  the  bed,  the  smaller  particles,  will  respond  slowly 
with  an  increase  in  response  speed  as  we  go  toward  the  bottom  of 
the  bed.   The  initial  response  of  bod  height  will  be  due  almost 
entirely  to  these  slowly  responding  particles  at  the  top  of  the 
bed,  but  there  will  be  considerable  particle  mixing  and  interac- 
tion within  the  bed  and  as  a  result  the  particle  distribution  at 
the  top  of  the  bed  will  be  continually  changing  and  yielding  dif- 
ferent response  characteristics.   It  must  be  pointed  out  that  this 
discussion  pertains  to  the  condition  at  the  top  of  the  bed  ^t 
times  before  t*.   After  t*  which  denotes  the  time  nt  which  a  dis- 
turbance is  to  re*ch  the  tor  of  the  bed  there  will  definitely  be 
a  change  in  roroslty  at  this  point  and  this  will  cause  varying 
response  characteristics.   When  considering  the  response  to  a 
step-down  the  situation  is  reversed.   The  particles  at  the  bottom 
will  still  respond  faster,  but  the  deviation  of  this  response  Is 
away  from  the  tor  of  the  bed.   The  particles  within  the  bed  will 
not  be  affected  greatly  by  particles  below  them  and  the  entire 
bed  wll]  settle  much  as  in  sedimentation. 

The  smaller  the  r^rtlcle  size  range,  the  more  neglible  the 
effects  described  above  would  be.   It  is  supposed  that  the  parti- 
cles used  by  Slis  et.  al.  (5)  were  very  uniform.   It  is  important 
to  be  conscious  of  the  Phenomena  described,  however,  since  In 
most  applications,  there  is  a  sizeable  particle  size  range. 

Effect  of  Aggregative  Fluidization 
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It  Mil  previously  mentioned  that  this  work  was  to  deal  chief- 
ly with  particulate  f luldization.      This  was   specified  since  the 
task  of  analyzing  the  unsteady   state  of  fluidlzation  processes 
would  he  greatly  magnified   if  non-idealities   were  present   in  the 
steady-st^te  condition  other  than  the  essential  non-linear  char- 
acteristics.     It   is   however  deemed   advisable  to  at   least   mention 
the  area   of  aggregate  fluidlzation  since  most   of  the  fluidlzation 
processes   fall    Into  this   category.      This   is  particularly  true  of 
gas-solid   fluidlzation  and    several  sources  have  recently  reported 
having  observ.d    it   in   liquid-solid   systems    (9).      Xa  liquid-solid 
fluidlzation  aggregative  fluidlzation   occurs  when  large,   heavy 
particles   are  fluidized  with  a  relatively  non-viscous   liquid  such 
as  water. 

In  this  present   investigation  a  tendency  toward  aggregate 
fluidlzation  was  observed   when   fluidizing  the  6   to   7  mesh  glass 
beads.      Two  distinct    ohases  were  observed.      One  consisted  of  areas 
densely  populated  with  particles   pnd   the   other  of    "bubbles"   of 
water  traveling  up  through   the-  bed.      Frequency  response   experiments 
were  carried   out  on  this   system   and  the  data   obtained  were  quite 
satisfactory.      A  problem  was   encountered   in  trying  to  apply  the 
linearized  -rodel    to   this   case,    howeve-.      The  terminal   falling 
velocity   of   the  particles  was  determined  and   n  calculated  by  the 
appropriate   correlation.      These  values  were  35.8  cm. /sec.   and    2.39 
respectively .      For  particulate  fluidlzation,    therefore,    one  would 
exoect  a   steady-state  relationship  from  Richardson  and   Zaki    (4) 
of  the  form 
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For  aggregate  fluidlzntlon  this  is  no  longer  valid.   Prom 
steady-state  data  obtained  for  this  sy^en  it  was  found  experi- 
mentally that  the  relationship 

4>S8  -  »-,*  ss  (65> 

could  be  used  to  represent  the  data.   For  this  case  a  was  found 
to  be  23.8  and  b  was  determined  as  2,$»      It  was  found  that  use 
of  the  values  of  Us  and  n  in  determining  the  system  time  constant 
gave  satisfactory  results  even  though  these  ere  not  the  same  as 
the  experimental  parameters  for  the  stendy-state  relationship. 
This  value  of  T  was  1.5  sec.   It  can  be  seen  from  Figure  (15)  that 
use  of  this  time  constant  gave  good  agreement  with  the  experimental 
data.   This  indicates  that  for  liquid-solid  aggregative  fluidized 
beds,  the  dilute  phase  does  not  contribute  significantly  to  the 
time  constant  of  the  system.   Although  there  is  no  theoretical 
development  to  surmort  this  procedure,  it  is  nrooosed.  as  a  possi- 
ble method  of  determining  the  transient  behavior  of  aggregative 
fluidized  beds. 

Application  of  the  Linearized  Model  to  Cther  Systems 

The  question  which  naturally  arises  at  this  point  is  whether 
or  not  the  theoretical  model  can  be  applied  to  systems  other  than 
those  on  which  experimental  d=ita  were  obtained.   This  is  obviously 
very  difficult  to  nredict  for  an  investigation  of  this  mture. 
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However,  If  It  is  realised  that  the  treatment  presented  is  an 
approximation  based  on  verified  steady-state  relationships,  it 
would  seem  th*t  the  model  could  be  expected  to  represent  the  re- 
sponse for  any  system  to  small  disturbances  where  reliable  steady- 
state  data  is  available.   It  is  porsible  that  the  model  could  be 
used  to  represent  gas-solid  systems  as  well  M  liquid-solid  sys- 
tems if  the  procedure  outlined  in  the  section  on  aggregate  fluid- 
ization  is  followed.   However,  much  more  experimental  work  is 
needed  before  the  general  applicability  of  the  model  is  proved. 

CONCLUSIONS 

The  following  conclusions  have  been  drawn  as  a  result  of 
this  investigation. 

1.  Fluidizatlon  processes  are  extremely  difficult  to  analyze 
from  a  dynamic  standpoint.   This  is  due  first  to  the  non- 
linear relationship  between  bed  porosity  and  fluidlzing 
velocity  and  second  to  the  different  mechanisms  for  bed 
extension  and  bed  contraction.   It  would  not  be  possible 
to  develop  a  single  model  to  represent  accurately  all 
transient  effects  of  a  fluldlzed  bed. 

2.  The  results  of  this  work  should  prove  useful  in  process 
control*   It  has  been  shown  that  the  response  of  the 
fluidized  bed  height  to  a  step  change  or  a  sinusoidal ly 
varying  fluidlzing  velocity  may  be  considered  as  a  first 
order  (time-constant)  system.   This  presupposes  of  course 
that  the  disturbance  encountered  will  be  small. 


7** 


3.   The  development  of  the  linearised  model  has  led  to  a 
theoretical  expression  for  the  time  constant  as  given 
by  Equation  (17)  or  (18). 

b.      Comparison  of  the  calculated  response  wit"  experimental 
results  show  that  the  limitation  of  the  linearized  model 
to  small  step  inputs  is  not  reflected  in  deviations  of 
the  calculated  time  constants.   It  is  seen,  however, 
that  the  size  of  input  affects  the  time  interval  over 
which  the  theoretical  model  is  applicable. 

5.  The  pulse  testing  method  of  obtaining-  frequency  response 
data  is  not  generally  useful  in  the  analysis  of  bed 
height  relationships  in  fluidized  beds.   It  is  recognized, 
however,  that  some  usefu"*  results  might  be  obtained  by 
using  a  more  efficient  pulse  form  than  the  rectangular 
pulse. 

6.  It  may  be  possible  to  treat  aggregative  fluidlzation  in 
the  same  manner  as  particulate  fluidlzation  in  determin- 
ing transient  bed  height  behavior.   It  is,  however, 
necessary  to  use  the  values  of  Us  and  n  to  calculate 
the  time  constant  for  this  tyre  of  system  even  though 
they  are-  not  the  experimental  parameters  of  the  steady- 
state  relationship.   This  makes  the  procedure  question- 
able. 

7.  Application  of  a  sinusoidal  input  to  a  fluidized  system 
will  cause  a  shift  in  the  steady-state  height  of  the  bed 
and  a  small  amount  of  distortion  can  be  expected  in  the 
output  wave  form. 
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8,  Particle  size  distribution  seems  to  have  a  definite 
effect  on  t he  dynamic  behavior  of  fluidized  beds. 
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NOMENCLATURE 

A     =  maximum  bed  expansion  (step  height),  cm. 

a     a  constant  in  Equation  (65),  cm. /sec. 

^col   ~  cross-sectional  area  of  column,  cm.  . 

^cy   ~  cross-sectional  area  of  cylinder  on  sine  wave  generator. 
cm.2. 

B  =  slope  of  linear  input  of  height  variable,  err. /sec. 

b  =  constant  in  Eauation  (65),  dimensionless. 

Dp  =  diameter  of  particle,  cm, 

h  =  bed  height  at  time,  t,  cm. 

hgg  -  steady-state  bed  height,  cm, 

h*  =  input  amplitude  in  bed  height  for  frequency  response,  cm, 

k  =  constant  in  Equation  (57),  sec,"  . 

L  =■  length  of  lever  arm  on  sine  wave  generator,  cm, 

n  -  constant  in  Equation  (1),  dimensionless. 

NSe  "  Reynolds  number,  dimensionless, 

rn  =•  height  of  staircase  function,  cm. 

t  =■  time,  sec, 

^0.02  =  time  interval  defined  by  Equation  (3*0,  sec, 

T  -  theoretical  time  constant,  sec, 

Us  a  terminal  falling  velocity  of  particles,  cm. /sec. 

V  «  volumetric  flow  rate,  ml. /sec. 

V^    =  tangential  velocity  of  lever  arm  on  sine  wave  generator, 
cm. /sec, 

Vs    -  velocity  of  oiston  in  sine  wave  generator,  cm. /sec. 
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Greek  Letters 

o(      -  phase  lag,  degrees. 

<-°   -  frequency,  1/sec. 

£   -  porosity,  dimenslonless. 

€«,_  s.  steady-state  porosity,  dimensionless. 

Gh  —  porosity  at  top  of  bed,  dimensionless. 

<J>   _  superficial  liquid  velocity,  cm. /sec. 

4>gs  -  steady-state  superficial  liquid  velocity,  cm. /sec. 

4?1  -   steady-state  superficial  liquid  velocity  after  step 
input,  cm. /pec. 

4>»  -  amplitude  of  sine  input  in  terms  of  superficial  velocity, 
cm. /sec. 


'  »  —     density  of  liquid. 
/° s  -  density  of  solid, 
•Y    =  viscosity  of  fluid. 
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APPENDIX 
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A.   Bibliographical  Essay 

At  the  present  time  the  need  for  experimental  and  theoretical 
studies  in  chemical  process  control  is  becoming  increasingly 
apparent.   In  the  past  ye^r  there  have  been  numerous  technical 
publications  on  the  response  of  chemical  reactor  systems  to 
changes  in  composition  and  temperature.   However,  the  pressure  and 
flow  fluctuation  and  their  interactions  with  composition  have  not 
been  treated.   Without  such  information  a  complete  analysis  for 
the  control  of  a  process  is  Impossible.   (10)   One  example  of 
such  a  process  is  the  response  of  the  liquid-solid  fluidized  bed 
to  changes  in  fluidlzlng  velocity. 

The  only  experimental  work  In  the  transient  behavior  of  fluid- 
ized beds  to  the  authors  knowledge  was  conducted  by  Slis,  et.  al. 
(5).   Their  work  involved  the  development  of  a  theoretical  ex- 
pression for  the  response  of  bed  height  to  a  step  input  in  fluid- 
izing  velocity.   Experimental  results  were  presented  to  verify 
their  equations. 

As  o^rt  of  a  study  on  the  feasibility  of  fluidized  bed  nuclear 
reactors,  Martin  Nuclear  Conroany  (8)  developed  a  differential  equa- 
tion to  represent  the  transient  behavior  of  bed  height  to  changes 
in  flow.   This  was  a  complicated,  non-linear  model  and  no  solu- 
tions were  obtained  in  analytical  form.   The  model  was  programmed 
on  the  analog  computer  and  responses  obtained  in  this  manner.   No 
experimental  work  was  conducted  in  their  investigation,  although 
the  authors  of  the  report  recognized  a  need  for  such  study  and 
so  stated  at  the  end  of  their  reoort. 
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It  is  evident  then  that  work  in  this  field  has  hardly  begun, 
although  there  is  a  definite  need  for  the  information.   Because 
of  the  interest  in  the  use  of  fluldlzed  beds  as  nuclear  reactors 
as  well  as  the  present  tendency  toward  more  accurate  process  con- 
trol, there  will  undoubtedly  be  more  investigations  of  the 
transient  behavior  of  fluldized  beds  in  the  ne?r  future. 
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B.   Procedure  for  Annllcation  of  Model 

The  following  shows  the  manner  in  which  the  curve  presented 
in  Figure  19  was  obtained.   The  problem  is  a  hypothetical  one 
which  deals  with  an  aspect  of  the  feasibility  of  fluidized  bed 
nuclear  reactors.   (8).   The  following  data  is  riven: 

€B8-  0.700 
h8g=  5.0  ft. 

cross-sectional   area   bed  -    3.^8  ft2, 
bed   diameter  -   7.0   ft. 
f  ss^    1.^0  ft. /sec. 
D     -    0.205   inches 
bed  weight  a    63,021  rounds 
temnerature   at  bottom  of  bed   -    472°F 
temperature  at  top  of  bed  =  574°F 
ave.    bed.    temperature   =   523°F 
f  f  ■=   48.2  lbs. /ft.3 
/0m'»   624.1   lbs. /ft. 3 
Us  -  3.29   ft. /sec. 
-4  f   -  0.025  cp.    -   1.68  x  10~5   lbs. /ft. -sec. 

NHe  ^    ^f^B    .    (48. 2)(0. 200(3.29)     _   1.61   x    i©5 
-7f  (12)(1.68)(l0-5) 

Since  NRe   is  well  over  500,    from   Equation    (32) 
n-  2.39 
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and 

the  steady-state  relationship  shou 
♦  as  =     use2s     =    3.29    e 

Id  be 

2.39 

ss 

T  is  calculated   from   Equation    (1?)   as 

y  ■=.           "ss 

(5) 

nOse^(l-C-se)        (2. 

39)(3.29)(0.700)1 

.39 

(0.300) 

-   3.48  sec. 

or 

from   Equation   (18) 

T   -                rSSeBS 

(CltO.MAl 

nC^ss(1-6ss)            (2.39)(1.4)(o.300) 

-  3.48   sec. 

Now 

the   steady-state  relationship  must 

be  linearized  about 

the 

point    (  <pBBt   €  as).      An  expression  for 

the   steady-state 

porosity 

is 

€       ^       i a_ 

SB                 x          TT 

hss 

whei 

-*»     a  —      bed  weight 

solid  density 

(63.000) 

(624.l)(38.48) 

cross-sectional  area   of  column 

=   2.63  ft. 

The 

steady-state  relationship  becomes 
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4>a.       ua(i.£i62)n 

'ss 


ss  h. 


differentiation  of  this  expression  gives 


12m  ,  ,itt.iA  "-1 

dh__  hsg 


^     n  U.    (l  -  M1)   n_1    U^D 


s 


which  Is  evaluated  at   h       -     5«0   to   give 

ss 


g     5S    =■    (2.39)    (3.29)  (1  -  *•&)!, 39  (Za&1\-  0.292   sec."1 
dhss  5.00'  '25.0' 


and  we  assume  the  straight   line  relationship 

hss         -  * 

-s&-      *        -     3.^3    sec. 

<p  ss  0.292 

The  first  input  to  be  considered  Is  a  linear  flow  input  of  4.62 

x  10"^   ftt  ,/sec.   Therefore  the  input  in  terms  of  bed  height  is 
sec. 

Mm   3.^3  <£  •=   (3.^3)(^.6?)(10"3)=  1.58  x  10-2  ft. /sec. 

and  the  input  Is  to  be  of  27  second  duration. 

The  differential  equation  giving  the  response  to  a  linear 
input  is 


dh  +   lh  ~  H 

dt      T     T 


which  hos  the  solution 
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-  | 
T 
h  -  h.fi  -     BT   (e      +     £     -  1) 

SS  m 


The  following   ooints   ^re  new  calculated   using  the  above  equation. 

t (sec. )  ss 

5  3.66  X  10-2 

10  1.06   X  10-} 

15  1.-3   X  10-| 

20  2.61  X  10~f- 

25  3.40  X  lo-J- 

?7  3.72   X  10-1 

At  this   point   the   inr>ut   is    to  remain   constant  at    (27   sec.)  (1.58) 

-2 
(10      )    ft. /sec.  -    0,426  ft,    for   10   seconds.      The  response  for  this 

portion   of   the  curve  w  ill  be 

-   (t-27) 

T  1 

h  -  Vi     =-  h  u(t-2?)   -   e  +  3.72   X  10"1 

s  s 


h=    0.426   -   0.372     -    0.044   ft. 

fcj     x  b-h    (ft.) 

t(sec. )  ss 


-2 


30  3.97  X  10 , 

35  4.12  X  10-} 

37  4.14  X  10-1 

Now  there  is  to  be  a  negative  linear  inout  with  B  "=■  1.58  X  10 

ft. /sec, 

.  (t-17) 

h  -  hss  =  BT  u(t-37)   e  -+-    {t"^i        -  1  +0.414 


This  int>ut  is  to  be  of  13.5  ^ec.  du^tion.   Calculated  values  are 
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h  -  h     (ft.) 
t(peot)  M 

ko  3.98  X  10"} 

45  3.36  x  10": 

50  2.62  X  10"* 

50-5  2.54  X  10"1 

Finally   the  input   is  to  be  constant   at   0.213   ft. 

..(t-50.5) 

h  -  h    =    (0.213-0.25^)   u(t-50-5)-e  4-      0.254 

s  s 

h  -  h     (ft.) 

4mJ  ss 

55  2.24  X  lO"1 

60  2.16  X  10-1 

65  2.14  X  10-1 

?0  2.13  X  10"1 

All  of  the  response  calculations   shown  here  are  plotted  on 
Figure  '9- 

C.       Outline  of  Proposed   Research  Project 

The  following  are  some  suggestions  concerning  further  re- 
search in  the  area  of  transient  behavior  of  fluidized  beds  and 
related   areas. 

1.      There   is   need  for  more  work  in  the  same  area  as   this 
present  work.      In  particular,    the  author  would  suggest 
that  additional   work  be  done  on  liquid-solid,   particulate 
fluidization  using  experimental  apparatus  with   a  greater 
range  than  the  apparatus   used  here.      If  a   sine-wave 
generator  were  constructed  so  that   the  frequency-range 
was  at   least   zero  to  two   radians  tier  second,    consider- 
able additional   data    could  be  obtained   on  phase   lag  and 
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amplitude  ratio  relationships  with  frequency.  Also  the 
use  of  moving  pictures  would  provide  an  extremely  accu- 
rate, but  relatively  expensive  method  of  analysis.  One 
of  the  better  methods  of  following  bed  height  is  present- 
ed by  Slis,  et.  al.  (5).  They  used  a  follow-up  system 
consisting  of  a  traveling  ohotoelectric  tube  and  a  ver- 
tical fluorescent  larr^. 

2.  The  dynamics  of  a  single  particle  suspended  in  a  flowing 
stream  is  possible.   To  accomplish  this  ■  column  at  least 
4  Inches  in  diameter  should  be  employed  for  the  study  of 
small  particles.   The  same  function  generating  elements 
that  were  used  in  the  fluldization  study  would  be  re- 
quired.  The  particle  could  be  suspended  by  using  a 
steady-state  flow  equal  to  its  terminal  settling  velocity. 
The  flow  could  be  changed  in  various  ways  and  the  motion 
of  the  particle  observed.   It  is  not  known  if  there  would 
be  any  oractical  value  for  such  work,  but  it  is  possible 
that  some  Information  concerning  basic  particle  dynamics 
could  be  gained. 

3.  A  study  of  the  effect  of  particle  distribution  on  the 
transient  behavior  of  fluidized  beds  is  recommended. 
Several  distributions  could  be  used  and  the  results  com- 
pared.  This  would  also  involve  a  fairly  extensive  study 
of  the  steady-state  axial  distribution  of  particle  dia- 
meter. 

**.   A  study  of  the  transient  behavior  of  semi-f luidized  beds 
would  be  useful.   Semi-fluidlzation  is  becoming  a  very 
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popular  nrocesp  nnd  there  will  surely  be  a  need  for 
work  in  this  field. 
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The  dynamic  response  of  bed  height  in  liquid-solid  fluidized 
beds  to  a  change  in  fluldizing  velocity  was  investigated.   This 
syster  is  recognized  as  a  very  non-ideal  system  from  the  stand- 
point of  dynamic  analysis.   Not  only  is  the  system  decidedly  non- 
linear, but  different  mechanisms  for  bed  expansion  and  contrac- 
tion are  evident. 

The  problem  was  at)r>roached  from  such  a  direction  as  to  lead 
to  information  useful  in  the  control  of  fluidized  processes. 
With  this  view  in  mind  a  linearized  model  was  developed  which 
could  be  expected  to  represent  the  response  of  fluidized  bed 
height  to  small  disturbances.   The  model  developed  was  linear  and 
first  order  expressed  as: 

dh    h    hi(t) 
dt    T      T 

where  h  is  the  instantaneous  bed  height,  h.  is  the  input  function 
in  terms  of  bed  height,  and  t  is  the  time.   Expressions  were  de- 
veloped such  that  the  time  constant,  T,  for  the  system  could  be 
estimated  using  the  parameters  characterizing  the  system: 

h  ss "ss  €  ss 


where  hss   is  the  steady-state  bed  height,   n   is  a   constant   dependent 
on  Reynolds  number,    U8   is  the  terminal  falling  velocity  of  the 
particles,   £  gs   is   the   steady-state  porosity,    and    & 8Q   is   the 


steady-state   superficial   fluid  velocity. 

Experimental  work  was   conducted   to  verify  the  linearized  mod- 
el.     The   inputs   considered  were  step,   pulse,   and   sinusoidal   in- 
puts  in  fluidizinp'  velocity.      The  agreement   of  the   linearized 
model  with  the  experimental  results  was  generally  quite  good  for 
each  of  these  inputs. 

The  response   to   a  step-up   in  fluidlzing  velocity  was   repre- 
sented much  more  accurately  by  the  linearized   model   than  the  re- 
sponse to  a   step-down.      This  was  because  of   the   tendency  of  the 
bed   height  to   give  a   linear  response  with  time  in   the  case  of  a 
step-down.      Based   on  these   considerations   the  time  Interval  over 
which  the  linearized  model  could  be  expected   to   give  an  arbitrary 
degree  of  accuracy  was   defined. 

An  attempt  was  made  to   obtain   frequency  resonse  data  from 
the  resDonse  to  a   pulse   Input  by  an  approximating  Fourier  trans- 
form method.      It  was   found  that  this  technique  could  not  be  de- 
pended  on  to  give  meaningful    results  when  applied  in  this   inves- 
tigation. 

The   results   of  frequency  response  analysis   showed  that  the 
distortion  appearing  in  the  output  wave  due  to  non-linearity  of 
the  system  was  not   severe  for  small  input   amplitudes.      In  this 
same  analysis  a   shift   In  the    steady-state  bed  height   was   observed 
when  a  sinusoidal    input  was  applied. 

It   is  proposed  that  particle  size  distribution  may  have  a 
significant   effect  on  the  transient   behavior  of  fluid ization  pro- 
cesses. 


It  was  found  that  It  la  possible  to  tre^t  the  dynamic  response 
of  aggregative  fluidlzed  bedB  In  much  the  same  manner  as  partic- 
ulate fluldlzatlon  was  treated  In  this  work. 

Transport  and  rate  processes  which  are  commonly  carried  out 
In  fluidlzed  beds  are  related  to  the  porosity  of  the  bed.  There- 
fore, no  study  of  such  a  orooess  would  be  complete  without  taking 
into  account  the  dynamic  behavior  of  porosity  or  bed  height. 


